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""^ ' A new, field-theory-based framework for discussing and interpreting exper- 

O" 1 imental tests of relativistic gravity, notably at the second post-Newtonian 

*h . (2PN) level, is introduced. Contrary to previous frameworks which attempted 

at parametrizing any conceivable phenomenological deviation from general 
relativity, we focus on the most general class of gravity models of the type 
suggested by unified theories: namely models in which gravity is mediated 
by a tensor field together with one or several scalar fields. The 2PN ap- 
proximation of these "tensor-multi-scalar" theories is obtained thanks to a 
diagrammatic expansion which allows us to compute the Lagrangian describ- 
ing the motion of iV bodies. In contrast with previous studies which had to 
introduce many phenomenological parameters, we find that the 2PN devia- 
tions from general relativity can be fully described by introducing only two 
new 2PN parameters, e and £, beyond the usual (Eddington) 1PN parameters 
/3 = (3 — 1 and 7 = 7 — 1. It follows from the basic tenets of field theory (no- 
tably the absence of negative-energy excitations), that /?, e and £ (as well as 
all the further parameters entering higher post-Newtonian orders) must tend 
to zero with 7. It is also found that e and £ do not enter the 2PN equations of 
motion of light. Therefore, second-order light-deflection or time-delay exper- 
iments cannot probe any theoretically motivated 2PN deviation from general 
relativity. On the other hand, these experiments can give a clean access to 7, 
which is of greatest significance as it measures the basic coupling strength of 
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matter to the scalar fields. Because of the importance of self-gravity effects 
in neutron stars, binary-pulsar experiments are found to constitute a unique 
testing ground for the 2PN structure of relativistic gravity. A simplified anal- 
ysis of current data on four binary-pulsar systems already leads to significant 
constraints on the two 2PN parameters: |e| < 7 x 1CT 2 , |£| < 6 x 1CT 3 . 

PACS numbers: 04.25.Nx, 04.50.+h, 04.80.Cc 
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I. INTRODUCTION 



The last three decades have been the golden era of experimental gravity: from Pound 
and Rebka to Hulse and Taylor, many complementary aspects of general relativity have 
been successfully tested. In particular, solar-system experiments allowed one to map out 
fairly completely weak-field gravity at the first post-Newtonian (1PN) approximation, i.e., 
to put stringent numerical constraints on a large class of possible deviations from gen- 
eral relativity at order 1/c 2 . Let us recall the useful role played in this respect by the 
first-order parametrized post-Newtonian (PPN) formalism |l|-[8| which introduced, in its ex- 
tended versions, about 10 independent phenomenological parameters to describe possible 
non-Einsteinian 1PN effects. Improved experiments are now planned to reach the second 
post-Newtonian (2PN) level -order 1/c 4 -, such as microsecond level light deflection exper- 
iments. Let us also mention that a 2PN treatment of the periastron advance is already 
significant for the binary pulsar PSR 1913+16 ||. It is therefore timely to undertake a 
systematic theoretical study of gravitational theories at this approximation. 

The 2PN limit of general relativity has already been studied in depth |T0|-[l8|,|9|,p!9| , but we 
also need to know what can be the possible deviations from these results in alternative the- 
ories of gravity. An ambitious program developed by Nordtvedt and Benacquista in [0-f22 



tries to extend directly the PPN formalism at order 1/c 4 , i.e., it aims at introducing a large 
number of parameters describing any possible relativistic theory at this order. Although it 
has only been partially implemented at the present time, this approach allowed one to derive 
some relations between these 2PN parameters by imposing the concept of "extended Lorentz 
invariance" (i.e., by requiring that the gravitational physics of subsystems, influenced by 
external masses, exhibit Lorentz-invariance) p^j22| . In spite of its partial achievements, 
the ability of such a general phenomenological approach to delineate the physically most 
important structures at the 2PN level is unclear. For instance, it was claimed in [^1[] that 
10 "parameters" are required to map Lorentz-invariant theories of gravity at the 2PN level; 
however, a careful reading of this article shows that several of these "parameters" are in fact 
functions of the distances between massive bodies, and could depend a priori on an infinite 
number of real parameters. 

In the present paper, we shall follow an entirely different methodology by developing 
a "theory-dependent" approach initiated in pM: Instead of considering any conceivable 
phenomenological deviation from general relativity, we focus on the simplest and best moti- 
vated class of non-Einsteinian theories, in which gravity is mediated by a tensor field (g^ u ) 
together with one or several scalar fields (y? a ). These "tensor-multi-scalar" theories arise 
naturally in theoretical attempts at quantizing gravity or at unifying it with other inter- 
actions (Kaluza-Klein and superstrings theories). Moreover, they are the only consistent 
field theories, containing only fields of infinite range, able to satisfy the weak equivalence 
principle (universality of free-fall of laboratory-size objects) ^[Q- in deed, massless gravita- 
tional theories incorporating, besides the metric g*^, vector fields, a second symmetric tensor 



Note, however, that the scalar couplings coming out naturally from unifying theories violate the 
equivalence principle |2G[ |. 
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field or an antisymmetric tensor field are known to present in general many flaws, such as 
discontinuities in the field degrees of freedom, negative-energy modes, causality violations, 
ill-posedness of the Cauchy problem, etc., not to mention the lack of theoretical motiva- 
tions for considering equivalence-principle-preserving couplings for such fields. By contrast, 
tensor-multi- scalar theories are well motivated, consistent and simple enough to allow their 



observational predictions to be fully worked out [^J]. Moreover, we believe that these field 
theories are the only ones satisfying the extended Lorentz invariance required in [^-|23|. It 
would be an interesting program to prove it rigorously. In that case, our field-theoretical 
approach gives a much more complete control of their structure than that of Refs. [2Tj-|2"3|], 
as exemplified by our results below. 

A detailed study of the 1PN limit of tensor-multi-scalar theories has been performed in 
24| , as well as the generalization of this approximation to the case of compact bodies (like 



neutron stars), called the first post-Keplerian (1PK) limit. We recall some of our results in 
section II below. Out of the 10 post-Newtonian parameters describing conceivable deviations 
from general relativity at the 1PN level, only two do not vanish in this class of theories: 
the parameters (3 = (3 — 1 and j = 7 — 1, introduced long ago (on different grounds) 
by Eddington |l|] 2 . In | 27| , j2"6 |, it was shown that the cosmological evolution generically 
drives these parameters towards values ^ 10~ 7 at our present epoch. This class of theories 
gives therefore a natural explanation (requiring no fine tuning nor the a priori presence 



of small parameters) to the bounds |t~| < 2 x 10 -3 [28] and \/3\ < 6 x 10~ 4 [29| found 
at the 1PN level in solar-system experiments, and furnish a motivation for increasing the 
precision of these measurements to the 10~ 7 level. Such an increase in accuracy down to a 
level comparable to 2PN effects (Gitiq/RqC 2 ~ 10 -6 ) makes it necessary to determine how 
2PN effects can influence such high-precision measurements, i.e., whether there are new 
and a priori unknown 2PN parameters which could complicate the interpretation of 1PN 
experiments. An example of this is given by higher-order light-deflection experiments which 



have been claimed to involve a new 2PN parameter [|I2],|T^|,|T3J . We shall, however, prove 
below that this claim is incorrect in the framework of tensor-scalar theories. 

The questions we shall address are thus: What are the new degrees of freedom describing 
the possible deviations of tensor-multi-scalar theories from general relativity at the 2PN 
order, and can the corresponding effects be separated from those associated with (3 and 7 ? 
On the other hand, do experimental bounds on j3 and 7 give constraints on possible 2PN 
non-Einsteinian effects ? 

Before entering into a detailed study of the 2PN limit of tensor-scalar theories, let us 
quote one of our main results: Two, and only two, new parameters arise at the 2PN level. 
We have denoted them by e and (, cf. Eqs. ( 3.3U| ) below. The possible 2PN deviations from 



the general relativistic physical metric tensor are given by 

,Ga{x!)U*W , 2C f.,G ff (x') t 3 , / G ( r(x")C/(x") 



2 The intuitively preferred role played by j3 and 7 in the PPN formalism is a further argument for 
working in the framework of tensor-scalar theories. 
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4 !( x)/,v^ + og,I) + o(I), 
*^(*) = o(£ ^)+o(±) , (Lib) 

^) = °(|'5) +0 (?) > ^ Llc ) 

where cx(x) is the mass density and U(x) = J d 3 x.' Ga(x')/\x — x'| the Newtonian potential. 
To increase the readability of Eqs. (1.1), we have suppressed the symbol ~ which should 
decorate all the quantities appearing in it: goo(x), U, G,a, . . .; see below. The same param- 
eters e, ( will be found to define the 2PN renormalizations of various Newtonian or 1PN 
quantities under the influence of self-gravity or external gravitational fields. 

In section II, we recall the action and the equations of motion of tensor-multi-scalar 
theories, as well as a few useful results concerning their 1PN limit. We also recall how the 
motion of self-gravitating bodies can be described in these theories. The main discussion of 
our paper, in section III, is devoted to the Lagrangian describing the motion of N massive 
bodies at the 2PN level. Our main technical tool is a diagrammatic expansion, which allows 
us to compute straightforwardly all the 2PN effects. In section IV we derive the 2PN metric 
corresponding to the Lagrangian of section III, and we verify and complement our results by 
considering the metric generated by one static and spherically symmetric body, whose exact 
solution has been derived in ||24j| . In section V, we discuss the impact of our findings on 



future relativistic experiments. We summarize our results and give our conclusions in section 
VI. To relieve the tedium, technical details are relegated to various appendices. Appendix A 
gives the explicit diagrammatic calculation of the 2PN Lagrangian. In Appendix B, we 
discuss the renormalizations of the Newtonian and 1PN coupling parameters due to 2PN 
effects. Finally, Appendix C derives the explicit 2PN formulae for the deflection of light and 
the perihelion shift of test masses. 



II. TENSOR MULTI-SCALAR THEORIES 

In this section we define our notation for dealing with tensor-scalar theories, and recall 



the results of [ 24] that we need below to study their 2PN approximation. 



A. Action and field equations 

For simplicity, we consider in the present paper only theories respecting exactly the weak 
equivalence principle, i.e., theories in which matter is universally coupled to one second 
rank symmetric tensor, say (^(a^). The action describing matter can then be written as a 
functional 

S m [lp m ,g^] , (2.1) 

where i[) m denotes globally all matter fields, including gauge bosons. Actually, from the 
perspective of modern unified theories, this class of models seems rather ad hoc. For instance, 
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string theory does suggest the possibility that there exist long-range scalar fields contributing 
to the interaction between macroscopic bodies, but all such scalar fields have composition- 
dependent couplings. However, a recent study of a large class of superstring-inspired tensor- 



scalar models ]26| has found that (because of deep physical facts) the composition-dependent 
effects represent only fractionally small (~ 1CT 5 ) corrections to standard post-Newtonian 
effects. 

At a fundamental level, the matter action S m should be chosen as the curved-spacetime 
version of the action of the Standard Model of electroweak and strong interactions, obtained 
by replacing the flat metric f^ u = diag(— 1, 1, 1, 1) by g^ v and partial derivatives by g- 
covariant ones. At a phenomenological level, the action describing a system of N (non 
self-gravitating) pointlike particles is 



N 

^2 / rn A cds A , (2.2) 

A=l ^ 



where ds A = [— ^^{x^dx^dx^] 1 ' 2 , and the m^'s denote the (constant) inertial masses of 
the different particles. The universal coupling to g^ implies in particular that laboratory 
rods and clocks measure this metric, which will therefore be called the "physical metric" 
[the names Jordan, Fierz, or Pauli metric are also used in the literature]. 

The difference with general relativity lies in that the physical metric 'g /iv , instead of 
being a pure spin-2 field, is in tensor-scalar theories a mixing of spin-2 and spin-0 degrees 
of freedom. More precisely, it can be written as 

g, v = A 2 (ip a ) 9 ; v , (2.3) 

where A(ip a ) is a function of n scalar fields (we choose A > to simplify some equations 
below). The dynamics of the pure spin-2 field g* u , usually called the "Einstein metric", is 
described by the Einstein-Hilbert action 

c 4 f d A x R* 
S sp in2 = t-tt —V9~*-Ti (2-4) 



AnG* J c 

where is a constant (the bare gravitational constant), R* is the scalar curvature of g*^ v 
(with the sign conventions of @]), and g* = \detg* \. On the other hand, the action 
describing the n scalar fields ip a reads 



c 4 f d 4 x ( 1 \ 
5spin0 = 4^T / ~^{~2 9 * Ulab ^ d ^ (padu(fb ) 



(2.5) 



where g% v is the inverse of g* the indices a,b,c, . . . vary from 1 to n, and 7 a ft(y2 c ) is a n- 
dimensional (a-model) metric in the internal scalar space spanned by the y? a, s. [7 a b must be 
positive-definite to get positive kinetic-energy terms.] A tensor-multi- scalar theory contains 
in general 1 + n{n — l)/2 arbitrary functions of n variables: the "coupling function" A(ip a ) 
involved in Eq. ( |2.3| ), and the n(n + l)/2 components of 7 a b from which must be subtracted 
n arbitrary functions parametrizing arbitrary changes of scalar-field variables ip' a = f a {(p b )- 
In the simplest case where there is only one scalar field (n = 1), the only arbitrary function 
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in the problem is A(ip), the unique component of the cr-model metric being always reducible 
to the trivial form 7 11 (y? 1 ) = 1. The reader should note that the consideration of multiple 
scalar fields, far from complicating uselessly our analysis, is in fact a technically powerful 
tool for delineating the structure of the possible deviations from general relativity. Once 
one is used to some notation, working with several fields is anyway not more difficult than 
working with only one. 

We could also have added a potential term V((p a ) in the action ( |2.5| ), but we will restrict 
our attention to infinite-range fields in the present paper (see |24[] for more details). Note 
that and the n scalar fields (p a are considered as forming the gravitational sector of the 
theory, by contrast with the matter sector described by the fields ip m of Eq. ( |2.1| ). 

Although one should always keep in mind that the metric measured by normal physical 
standards is g^, it will be convenient in the following to formulate the theory in terms of 
the pure spin-2 and spin-0 fields g* u and ip a . The Einstein- frame infinitesimal lengths £, 
time-intervals t, and masses m will therefore be related to the physical (measured) ones by 

t = A~\tp)7 , t = A' 1 ^)!, m = A{ip)m . (2.6) 

For instance, the action ( |2.2|) describing N (non-self-gravitating) pointlike particles can be 
rewritten as 

N . N . 

S m = - ^2 / mAcA((p a (x A )) ^-g* l/ (x A )dx l ^dx A = - ^ / m A (ip a (x A )) cds* A , (2.7) 

A=l J A=l J 

where the Einstein-frame masses m A ((p a ) = A(ip a )m A are no longer constant, as opposed to 
the m A s. 

The field equations deriving from the total action S sp i n 2 + Sspino + S m read 

r; u = 2 lab ( v ) d^d vV h + ^ (t; v - l - t* 9 ;)j , (2.8a) 

a {g ^ a = -^a a (p)T* , (2.8b) 

SS m [lpm,gixu]/Slp m = • (2.8c) 

Here R*^ is the Ricci tensor of g*^, = (2c/ y^) 5S rn [ip m , A 2 g*^)/6g*^ is the Einstein- 
frame energy tensor (related to the conserved "physical" energy tensor T^ LU by T£ v = A Q T^ U , 
see plf), and the d'Alembertian □( fl *, 7 ) is covariant with respect to both space-time and cr- 
model indices, i.e., involves the Levi-Civita connections of both g* and •y ab (ip), denoted 
respectively as T*J^ and 7fe c (y?) : 

U {g ^ = £> [d,d^ a - V; x u d xV a + ^Md^d^] . (2.9) 

In Eqs. ( p.8|) and everywhere else in this paper, the various indices are moved by their 
corresponding metric, for instance T* v = g* a gtpT^, T* = g^ v T^ v , = g m T au , etc., but 
also a a (ip) = -y ab ((p)a.b(tp) where ^ ab is the inverse of 7 a b and where we have introduced the 
notation 
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(2.10) 



Note that, in view of the third equation fl2.6|) , the definition of ct a (ip) can be rewritten as 

aM = , (2.11) 

where m(<p) is the mass of any non self-gravitating particle in the Einstein conformal frame. 
This second way of defining a a is quite general as it encompasses both self-gravity effects 



(see below) and possible composition-dependent effects [f26 |. 



B. 1PN approximation 

As is clear from Eq. ( |2.8b| ), the quantity a a (ip) (which is a vector field in the internal "<r- 
model" space spanned by the (p a, s) plays the crucial role of measuring the coupling strength 
of the scalar fields to matter. As we shall see below, all post-Newtonian deviations from 
general relativity (of any post-Newtonian order) can be expressed in terms of the asymptotic 
value of a a (ip) at spatial infinity (i.e., far from all material sources) and of its successive 
scalar-field derivatives. Denoting by D a the covariant derivative with respect to the internal 
metric •y ab , we define 

dot 

(3 ab = D a a b = — - j c ab a c , (2.12a) 
f3' abc = D a D b a c . (2.12b) 

Denoting by v?o the (cosmologically imposed) background values of the scalar fields, we then 
set a° a = a a (ip ), (3 Q ab = (3 ab (^p ), (3® hc = (3' ahc (yo), etc., the index o always meaning that these 
cr-model tensors are calculated at tpo. As shown in [Q, the effects of the scalar fields on 



any observable effect at the first post-Newtonian level depend only on two contractions of 
a a and Pah, namely^ 

a 2 Q = a° a a a = a^^o)^ , (2.13a) 
(a(3a) = a a (3° ab a b . (2.13b) 

Indeed, the effective gravitational constant between two massive particles is given by 

G = G,A 2 (l + a 2 ) , (2.14) 



(with A = A((p )) instead of the bare constant G* involved in the action ( |2.4j ), and the 
Eddington parameters ft, 7 read 



3 Note that in the simplest one-scalar case tp = ip 1 , these contractions reduce to simple products, 
«o = a o x a 0i (otfta)o = o-q x Pq x ao, where oq = a(po) = dkiA((po)/d(po, and fto = da(ipo)/dp>o. 



8 



7 = -2aj|/(l + ajj) , (2.15a) 
P = ^(af3a) /(l + a 2 ) 2 , (2.15b) 

instead of their general relativistic values (3 = 7 = 0. Let us recall again that (3 and 7 
are usually denoted as (/3 — 1) and (7 — 1) in the literature. We will avoid the notation 
(3, 7 in the present paper to prevent a possible confusion with the cx-model tensors (3^ and 
■jab- Note, however, from Eqs. ( j2.13| ) and ( |2.15|) , that our notation has been chosen so that 



7 oc ^abCt a cc and (3 oc (3 a ba a a . As shown in [Q, the results ( |2.14| )-( |2"T5[ ) can be simply 



interpreted (and remembered) in terms of the exchange of gravitons and scalar particles 
between material sources. For instance, Eq. (|2.14j ) is the sum of the usual contribution of a 
graviton exchange, together with the contributions of scalar exchanges, X]a^* a a a 0' see 
Fig. |I|. The global factor Aq is due to the change of units between the Einstein metric and 
the physical one = A 2 ((p)g^ u used to measure forces (see Eqs. ( PTB| ) above). Similarly, the 
parameter (af3a)o involved in (|2.15b| ) corresponds to an exchange of scalar particles between 
three massive bodies, as shown in Fig. ^. The method that we will use in section III below 
to study the second post-Newtonian approximation is a straightforward generalization of 
these diagrammatic observations. 



C. Self-gravitating bodies 

When studying the motion of massive bodies in the solar system, several dimensionless 
ratios happen to be small. If we denote by m, v and R the typical mass, orbital velocity 
and radius of a body, and by r the typical distance between two bodies, we find Gm/rc 2 ~ 
v 2 /c 2 < 2 x 1(T 8 for the fastest planets, while Gm/Rc 2 ~ 2 x 1(T 6 for the Sun and ~ 7x 1(T 10 
for the Earth. This is the reason why the formal "post-Newtonian" expansion in powers of 
1/c 2 is so useful for analyzing the predictions of relativistic theories of gravity in the solar 
system. However, in situations involving compact bodies, like neutron stars (pulsars) in 
binary systems, it is necessary to distinguish the self-gravity parameter s ~ Gm/Rc 2 ~ 0.2 
from the orbital parameters Gm/rc 2 ~ v 2 /c 2 <C 1. In that case, one can still describe the 
motion of the bodies by means of an expansion in powers of Gm/rc 2 ~ v 2 /c 2 ~ l/ c orbitai> but 
one must not expand in powers of the compactnesses s. Such an expansion scheme has been 
called "post-Keplerian" in [p4j] , since it is closely linked with the phenomenological approach 



to the analysis of binary pulsar data, introduced in |3T,32|] under the name of "parametrized 
post-Keplerian" formalism. 

In the present paper, our main goal is to analyze tensor-multi-scalar theories of gravity 
at the second post-Newtonian level, i.e., including all terms of formal order 1/c 4 (be them 
of "orbital", "self-gravity" or mixed type). We found that the most efficient way of doing 
so is to derive first the second post-Keplerian (2PK) limit of these theories by means of a 
diagrammatic method. The 2PN approximation is then obtained by expanding our general 
2PK results in powers of the compactnesses, up to the required order in s. 

Let us recall how the motion of self-gravitating bodies is described in tensor-scalar 
theories. Following a suggestion of Eardley |53 |, one skeletonizes extended self-gravitating 



bodies as pointlike particles whose inertial masses m^ depend on the scalar fields (p a , as 
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opposed to the constant m^'s in Eq. (|2.2|) describing non-self-gravitating bodies. This scalar 
dependence of the inertial mass is due to the influence of the local scalar field background on 
the equilibrium configuration of the body. The action describing N self-gravitating bodies 
is thus written as 

N r N 

Sm = ~J2 ^A(^)cds A = - X) / m A^ a )cds\ , (2.16) 
A=l A=l J 



where m^(y) = A({p)m A ({p). The validity of the skeletonized action (|2.16|) has been justified 
in Appendix A of pi by a matching argument. 



Note that the second expression ( 2.16Q , in terms of the Einstein line element ds A , is for- 



mally identical to Eq. ( |2.7| ) describing non-self-gravitating bodies. However, the important 
difference is that m A (ip) can now be a non- universal (body-dependent) function of the scalar 
fields, instead of being merely proportional to A(<p). It is therefore convenient to generalize 
to the case of self-gravitating bodies the a-model tensors defined in Eqs. ( |2.10| )- (p.l2j) above: 



a _ dhi m A (<p) ( , dlnm A (y) 

«• = dv a = <<P) + > ( 2 - 17a ) 

(3i = D a at, (2-lTb) 
f3' a i = D a D b af , (2.17c) 

and similarly for higher covariant derivatives. (As above, we will raise and lower the a-model 
indices a,b, . . . with r y ab or 7 a &.) Here again a a A plays the fundamental role of measuring the 



coupling strength of the scalar fields to the self-gravitating body A. Indeed, equation (|2.8b|) 
now reads 

^/^-^E^, (2.18) 

° A 

where T A = —rriAC 2 (ds A /dx ) 5^ 3 ^(x — x^) / y g*(x A ) is localized on the position of the A- 
th "particle". Note that the body-dependent quantities ([2.17]) reduce to the definitions 



( |2.10] )- (|2.12| ) for non-self-gravitating bodies, since the inertial masses m A are constant in 
that case. 

As shown in [Q, the 1PK approximation of tensor-multi-scalar theories can then be 
expressed very simply in terms of contractions of and f3^ b . More precisely, the Lagrangian 
describing the motion of iV (spherical) self-gravitating bodies at the 1PK level reads (in 
Einstein-frame units) 

A A^B B^A^C ^ ' 

where the notation B 7^ A 7^ C excludes A = B and A = C but not B = C, and where 

L« = -m° A <?yfl-vy<* = -mY + \m A y\ + + O (1) , (2.20a) 
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T (2) _ G AB m A m° B r. . 3 , , 1 



~ 7"AB 



1 + 2? K + V ^ )- 2? (Va ' Vb) 
-^(n AB -v A )K B -v B ) + %(v A -v B ) 2 , (2.20b) 



2c 2 c 



c 2 

lTc = - G ^ G ^\ m y i + A 

tabTacc 1 

Here we have set m A = m^(ipo), and we have denoted by thab = yabI t ab (with tab = 
— xb) the unit vector directed from body B to body A. This Lagrangian involves three 
body-dependent parameters generalizing the effective gravitational constant ( |2.14 ) and the 
1PN Eddington parameters ( |2.15| ), namely 

G AB = G.[l + {a A a B ) } , (2.21a) 

7y = _ 2 J^L ( 2 21b) 

B 1 + {a A a B )o 
Pbc - - r 7 ^f c \ , (2.21c) 

whereQ (a^ag) = a A "y a ba b B , {olb^aolc) = a B^ab a c^ an< ^ the m dex o means that these con- 
tractions are calculated at <p a = (p^. Of course, when using physical units related to the 
asymptotic metric g^ u , the effective gravitational constant reads Gab = C*Aq[1 + (q;a«s)o]; 

A 

and the parameters 7^5, Pbc d° no ^ change since they are dimensionless. Note that Figures 
HI and ^|give again a diagrammatic interpretation of Gab and (qbPa<^c)o, with the coupling 
coefficients a a and (3 ab of each body being replaced by their strong-field counterparts a\ and 

Pa- 

The pointlike description ( |2.16|) neglects all finite size effects. One might a priori worry 



that scalar interactions might introduce relativistic couplings to "scalar" multipole moments, 
starting with the spherical inertia moments / ~ J (i 3 xcx(x)x 2 . If that were the case, that 
would introduce in the equations of motion fractional corrections of order 7(1; /c) 2 (R/r) 2 to 
the leading Newtonian force. Such terms would be important notably in the Moon-Earth 
relative motion. In fact, no such terms exist. This can be seen in two ways: either by exam- 
ining the exact external tensor-scalar field of a spherical extended body (see below) which is 
found to depend on the structure of the body only through tua and a a, or from considering 
the general multipole expansion of the retarded field generated by an arbitrary scalar source 
S(x) which contains only one monopole moment M (u) — ~ J d 3 x dzS(x., u + zr/c) and 
higher (i > 1) multipole moments (see e.g. ||34|| ). Because the bodies making up the solar 
system are nearly spherical and in inner equilibrium^, one easily checks that the corrections 



4 Again, in the simplest case of one scalar field ip = ip 1 , these quantities reduce simply to a a x ots 
and as x @a x «Ci where o<a = dlnrriA/dip, (3a = dctA/dp. 

5 However, if one wanted to deal with spherical bodies which are not in equilibrium, say a collapsing 
star, one should carefully take into account the effects associated with the intrinsic time variation 
of the scalar monopole moment m^a^. 
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induced by scalar-mediated couplings to the scalar multipole moments of order i > 1 are 
negligible compared to the interactions described by the action ( 2.16 ). As for the correc- 
tions induced by tensor-mediated couplings to the mass and spin multipole moments of order 
t > 1, we assume that they are properly added, following for instance the recent papers |35 
which worked out a general relativistic description of multipole interactions. 



III. iV-BODY LAGRANGIAN 



Before entering the technical details of our derivation of the iV-body Lagrangian, let us 
clarify the approximation methods we shall employ. As discussed in the previous section, 
we shall first study the second post-Keplerian approximation of tensor-multi-scalar theories, 
before particularizing our results to the second post-Newtonian case. In other words, the 
compactnesses s ~ Gm/Rc 2 of the bodies will not be considered a priori as small parame- 
ters. In order to construct the Lagrangian describing the motion of N self-gravitating bodies, 
we will eliminate the field degrees of freedom, i.e., we will solve for the Einstein metric g*^ 
and the scalar fields ip a in terms of the material sources, using the field equations ( ^.8a| ), 
( p.8b|) . To perform this elimination, we shall consider that the interaction is propagated by 
a time-symmetric (half-retarded-half-advanced) Green's function. This leaves out radiation 
damping effects. However, the latter are negligible when studying weakly self-gravitating 
bodies at 2PN order. Indeed, in general relativity, the leading dissipative effects occur 
at order 1/c 5 , because of the well-known quadrupolar radiation of gravitational waves. In 
tensor-scalar theories of gravity, the leading radiation emitted by systems of compact bodies 
is dipolar and occur at order s 2 /c 3 . In the solar system case, s ~ 1/c 2 and the dipolar radi- 
ation is negligible compared to the 0(l/c 5 ) effects due to the spin-2 quadrupolar waves and 
spin-0 monopolar and quadrupolar waves |[24|| . In either case, we compute the conservative 
part of the gravitational interaction and assume that (tensor and scalar) radiation damping 
effects are added separately when they are not negligibly small. 



A. Diagrammatic expansion 



We want to construct a "Fokker" Lagrangian |36| , p!7| describing the motion of N massive 
bodies by eliminating the field degrees of freedom from the total Lagrangian of the theory. 
In order to prove formally that such a construction reproduces the correct dynamics of 
the bodies, let us introduce a global notation for the fields <3> = (g* — f^ u , (p a — where 

= diag(— 1, 1, 1, 1) is the flat metric and y?Q are the background values of the scalar fields. 
Similarly, we denote globally by a the matter variables, i.e., the N massive worldlines xa 
involved in the matter action ( [2.16j ). The total action of the theory can therefore be written 
in the schematic form 



S tot [*,<S>] = S$[<S>) + S m [a,<S>) , (3.1) 

from which we want to eliminate the fields $ by expressing them in terms of the matter 
variables a. This elimination cannot be done when working directly with the Einstein- 
Hilbert action (|2.4|) because of its invariance under diffeomorphisms (technically, the kinetic 
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term of the gravitons is non-invertible). We need to reduce the field equations by means of 
a specific coordinate condition in order to solve for g* . In the language of particle physics, 
we need to fix the gauge in order to define the propagator of the gravitons. We will choose 
the g*-harmonic gauge for our explicit calculations. We then replace by its gauge- 

fixed version S% '[<&] = S$ [$] + gauge- fixing terms (see Appendix A below), and the field 
equations deriving from the total action S^t [°~-> read: 



6S*£[v,$] _ 5 Sf ■{<!>] | 6S m [<T,*] 
5$ 5$ <5$ 

;.f. 



, (3.2a) 



SJ&S = SSrn^S = o . (3.2b) 
5a 5a 

Eq. (|3.2a|) can now be solved perturbatively, and we denote by $[cx] its solution. The Fokker 
action (which is a functional of the matter variables only) is then defined as 

S F [a] = [a,<§[a]] = Sf [%]] + S m [a,$[a\] , (3.3) 

and its variation with respect to a reads 

5S F [a] (8S&[v,*]\ , (5S&'[a,Q]\ x 6*\a\ _ (g _ 4) 




5a \ 5a I _ \ 5<f> I _ 5a 

$ = $[cr] 

Since $[er] is precisely the solution which annuls the functional derivative 5Sf ot '/5&, the 
second term on the right-hand side vanishes, and one finally gets 



5S F [a] (5S^[a,^}\ _ ( 5S m [a,$\ 



5a \ 5a I _ \ 5a / .., TA , 



(3.5) 



Comparing ( p.5|) with (|3.2b|) , we see that we have formally proved that the Fokker action 
( |3. 3Q gives indeed the correct equations of motion, i.e., describes the actual dynamics of the 
matter variables a in presence of the fields $[<t]. It should be noted that S F [a] is not simply 
given by the matter action S m [a,<&[a]\ embedded in the self-consistent background field 
but that the field action S$ L \&[a]\ also contributes to the dynamics of the material 
bodies. 

To start with, we do not need to assume that the bodies are moving slowly with re- 
spect to the velocity of light. In technical terms, we shall expand the Lagrangian in powers 
of G, i.e. Gm/rc 2 ( "nonlinearity expansion"), while keeping unrestricted the magnitudes 
of v/c and s. We need to retain the terms up to order G 3 in the nonlinearity expansion 
of the Fokker action ( |3.3| ). The zeroth order term in S F is of course the action describ- 
ing free bodies So[a) = S m [a, $ = 0], i.e., the matter action computed for g^ u = 

and (p a = (Pq. It can be written explicitly as So = YIaI dtL^, where was given in 
Eq. (|2.20a|) and has the structure mc 2 (l + v 2 /c 2 + t> 4 /c 4 + t> 6 /c 6 + •••). The next approx- 
imation, first power of G, has the structure mc 2 x (Gm/rc 2 ) x (1 + v 2 /c 2 + t> 4 /c 4 + • • •) 
and describes the two-body interaction, starting with the Newtonian term G AB m A m B I t ab 
(with the self-gravity- modified effective gravitational constant Gab of Eq. ( |2.21a| )). The 
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second power of G, mc 2 x (Gm/rc 2 ) 2 x (1 + v 2 /c 2 + •••), describes three-body interactions 
starting with the self-gravity- modified 1PK term oc G abG Ac m A m B m c I r ab t acc 2 ■ Finally, 
the G 3 level corresponds to four-body interactions: G 3 mAfriBTnci r n£,/r 3 c 4l + self-gravity and 
velocity modifications. Instead of counting powers of G, we can count the number of matter 
source terms involved: we must keep up to four powers of a, i.e., of the masses. Since the 
matter action S m [a, <£>] is linear in a (see Eq. ( |2.16|) ), the solution $[er] of the field equations 
( |3.2a| ) starts at order a. We need therefore to expand the total action S^ t '[a, $] up to orders 
0($ 4 ) and 0(cr$ 3 ) included, before replacing $ by its solution $[cr]. 
Let us first expand Sf ot ' in powers of the fields $, and define 

S^ t [a, $] = S [a] + St[a, $] + S 2 [a, $] + S 3 [a, $]+••-, (3.6) 

where the term Si[a, $] involves the i-th power of <3> (and zero or one power of the material 
sources a). For instance, Si is the linear interaction term between the fields and matter, 
and has the formal structure S\ = aa$>, where a is a coupling constant. On the other 
hand, the term quadratic in $ has the form S 2 = — + i/3<r$ 2 , and involves both the 
kinetic operator oc □ (or "inverse propagator" V^ 1 ) of the fields and a vertex describing the 
interaction of matter with two fields (with a coupling constant 0). It will be convenient to 
introduce a diagrammatic notation for this expansion. [Bertotti and Plebanski ]37j were the 
first to introduce a similar diagrammatic notation for solving Einstein's equations perturba- 
tively.] Let us denote the propagator V by a straight line, the material source a by a white 
blob, and the term (T 3-1 "^) by a black one; see Fig. || As this diagrammatic representation 
will be an important tool in the present paper, let us explain it in detail with a simpler 
example. 

Let us consider the action (in Minkowski spacetime) 



S[ip] = / d 4 x 



(3.7) 



where (dtp) 2 = f^d^ipdvip and where a(x) is a given (spacetime distributed) source for <p(x). 
Integrating by parts, the kinetic terms for <p(x) read +(l/87r) f d A x tp(x)Of(p(x), where □/ 
is the flat d'Alembertian. Identifying this with — ^ipV~ 1 f [which is a symbolic notationF] for 
— I ff d 4 x d 4 y ip(x)V~yip(y), where V~y is the kernel of an operator acting on functions ofx M ], 
we get T~ y x = — (47r) _1 n :!; (5^(x— y). The inverse of the operator V~ y l is the propagator V xy = 
Q(x — y), where Q(x — y) is a (translation-invariant) Green function, solution of U x Q(x — y) = 
— Atc5^(x — y). The cubic vertex V3 is defined as the distributional kernel entering the term 
3S3 = / dx g(p(dip) 2 , i.e., 3S3 = j dx\ dx^dx^V^Xx, X2, Xs)<f(xi)<f(x2)<f(xs). Requiring this 
kernel to be symmetric in its arguments leads to the explicit expression 

g r d d d d 

V 3 (X 1 ,X 2 ,X 3 ) = - [g^I S ( X l ~ X l)-Q^Jlb( X l ~ X 3.) + Q^ S ( X2 ~ X 3)g^jI S ( X 2 ~ X l) 

d d 1 

5(x 3 - xx)— T 5(x 3 - x 2 ) , (3.8) 



dreg v 17 dx^ 



6 Note that in the operator notation used here, any "contraction of spacetime indices" means an 
integration over the corresponding spacetime coordinates, e.g. (V<p)(x) = f d A yV xy tp(y). 
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where the factor 1/3 comes from the average over the three different permutations 
needed to symmetrize V 3 (xi,x 2 ,x 3 ). Similarly, the quartic vertex, defined by 4S4 = 
f dxidx2dx 3 dx4Yi(xi, x 2 , x 3 , X4)(p(xi)ip(x2)(p(x 3 )(p(x4) , is 

V i (xi 1 x 2 , x 3 , x 4 ) = X5(xi - x 2 )5(xi - x 3 )5(xi - x A ) . (3.9) 

In the diagrammatic notation of Fig. ^|, the blobs denote some spacetime functions (cr(x) 
for the white blob and — (4ir)~ 1 0(p{x) for the black one), and a line denotes a propagator 
*Pxy = G( x ~ y)- Connecting a line to a blob or to a vertex (which is a cluster of several 
infinitesimally close points) means that one "contracts" {i.e., integrates) over the points at 
the extremities of the line. For instance the T-shaped diagram on the left of the third line 
of Fig. [7| below would represent, in the model ( |3.7|) . 



dx 1 dx 2 dx 3 dy 1 dy 2 dy 3 V 3 (x ll x 2 , x 3 )g(x! - yi)G(x 2 - y 2 )G(x 3 - y 3 )cr(y 1 )a(y 2 )a(y 3 ) 

f 3d 

= g dxdy 1 dy 2 dy 3 G{x -ytj—Gix -y 2 )-^G(x -y 3 )a(y 1 )a(y2)(T(y 3 ) . (3.10) 

Similarly, the X diagram on the left of the last line of Fig. [7] would denote 

dx 1 dx 2 dx 3 dx 4 dy 1 dy 2 dy 3 dy 4 V 4 (x 1 , x 2 , x 3 , x 4 )Q(x 1 - yx)Q{x 2 - y 2 )Q(x 3 - y 3 )G(xi - y 4 ) 
x ^lM^M^M^) 

A / dxdyidy 2 dy 3 dyiQ(x - yi)Q{x - y 2 )Q{x - y 3 )Q{x - y i )a{yi)a{y 2 )a{y 3 )a(y i ) . (3.11) 



Having explained the precise meaning of our symbolic notation, let us come back to the 
general action (|3.6|). 

The different terms of the expansion ( |3.6| ) can be represented as in Fig. |], which defines 
the different diagrams. Note that, as in our example, we have conventionally factorized a 
coefficient 1/i in defining the i-linear vertex Vi from the 0($ l ) action: Si = V^ji. This 
coefficient is chosen in order to simplify the field equations (|3.2a ), whose diagrammatic 



expansion is displayed in Fig. |5|. [The reader is invited to derive for himself the field equations 
of Fig. [5] from the action of Fig. ^, keeping in mind that the multilinear forms in $ appearing 
in Fig. I are supposed to be symmetric in <&(xi), . . . , Thanks to Euler's theorem on 

homogeneous functions, the field equations imply the useful result 

$ x 6S ' ot }^ $] = Sx + 2S 2 + 3S 3 + ■ ■ ■ + iSi + ■ ■ ■ = . (3.12) 

In diagrammatic terms, this corresponds to inserting a black blob at the free ends of the prop- 
agators in Fig. |, i.e., to express the kinetic term of the fields, ^&V~ 1 & = ^(V^^ViV^ 1 ^), 
in terms of the other diagrams. 

The Fokker action ( |3.3|) can now be written straightforwardly by replacing in the total 
action ( |3.6| ) the fields $ by their solution $[<r], i.e., by replacing the black blobs in terms of 
the white ones through an iterative use of Fig. |[ The most delicate term to compute would 
be the contribution due to the kinetic term of the fields in S 2 , because one must expand up 



15 



to order a 3 the two fields <£> it involves. Fortunately, one can avoid estimating this term by 
using the Euler identity ( p. 12 ) to eliminate it from the Fokker action: 



S F [a] 



(S + S 1 + S 2 + ---)--(S 1 



S + 



2 2 6 



2S, + 3S* 



0(cr z 



*=*[<t] 



(3.13) 



The result of inserting Fig. |5| into Eq. ( |3.13| ) is displayed in Fig. |6]. [The different diagrams 
have been drawn so that angles appear only at the vertices involving matter sources.] In the 
following, we will designate these diagrams by the letter they most naturally evoke, so that 
the final result for the Fokker action reads 



Spier] = SnW] 



Ii 

2 



iv + i T 

2 3 



3 2 



I H + ix 

2 4 



0(a 5 



(3.14) 



The explicit form of this action is now only a question of straightforward algebra: one must 
expand Sf ot ' [a, $] up to order <3> 4 included to get the expressions of the field propagator V 



and of the vertices defined by Fig. f|, and merely replace them in Fig. [6|, i.e., Eq. fl3.14| ). 
In doing so, one must take into account the fact that $ is a global notation for both the 
gravitons h^ u = g^ u — f^ u and the scalar fields ip a — </?o, and therefore that each propagator 
link in Fig. ^ is to be replaced by a sum of terms corresponding to the different fields. To 
simplify the expansion of the scalar-field action (|2.5| ), it is convenient to choose Riemann 
normal coordinates at ip^, so that the metric 7 a ft(v 9 ) can be written as 

labfa) = labifo) +0X(/- if*) - -Rad,d{<PoW ~ <fo) (/ ~ <fo) + 0{? ~ Vof , (3.15) 

where R a bcd is the Riemann curvature of 7^. This choice cancels the term of order tpdipdtp in 
•S'spinO) *- e -) the "T" vertex connecting three scalar fields. The different diagrams of Eq. (|3.14|) 
can thus be decomposed into the elementary diagrams displayed in Fig. 0, where curly and 
straight lines represent the graviton and scalar propagators respectively. The coefficients 
appearing in this figure are simple binomial coefficients coming from the various ways of 
choosing the lines (see Appendix A). The diagrams involving only gravitons give the 0(G 3 ) 
approximation of general relativity, which has been studied in the literature (see notably 



14,17]). The diagrams involving at least one scalar propagator give therefore all the looked 



for deviations from general relativity predicted by tensor-scalar theories at this order. Their 
explicit calculation is performed in Appendix A below. Their global structure is easy to 
grasp. Denoting by Q a Green function, solution of OfQ(x) = —4k6^(x) where □/ is the 
flat d'Alembertian, each scalar propagator is V® b oc Q{xa — XB)l ab , while each graviton 
propagator is V^ s oc Q(x A - x B ){f ai fps + fasffr - fapf-ys)- Eacn matter vertex containing 
a single scalar coupling brings a factor a^, those containing double scalar couplings bring a 
factor f3^ h + a^a^, while triple scalar couplings bring a factor (3'^ bc together with extra (3^ h 
and terms [see Eq. ( |A11| )-( JA14[ ) in Appendix A below]. Each link implies a contraction 
over the internal indices^. Finally, the global structure of the action is 



7 Those internal contractions make it very easy to work with n scalar fields. Actually, once one is 
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A B •* 

x /K B ,J< B Jvo) 1 ftjVo),C,...fe),i? aW W) , (3.16) 



where is the Minkowski proper time along the worldline x A (r A ), and u\ = dx A /cdr A . We 
use in our calculations the symmetric Green function |(£? re tarded + ^advanced), given by 



Gsym(x A -X B ) = 6 (fnv(x^ ~ X B )(x A - X B )) 

= [5(x° A -x° B -r AB ) + 5(x° A -x° B + tab)] /2r AB . (3.17) 

The power of our diagrammatic approach shows up in the fact that one can identify 

the theory-dependent parameters appearing in each term of the action without doing any 

calculation. For instance, the interaction between two bodies A and B is described by 

the first two "I" diagrams of Fig. |7|, and involves therefore (because of the first diagram) 

the contraction (a^agjo- This explains why the effective gravitational constant G AB and 

the 1PK Eddington parameter ^y AB appearing in L^ AB , Eq. ( |2.20b|) , depend only on this 

contraction. Similarly, the interaction between three bodies is described by the "V" and 

"T" diagrams of Fig. [7|, and we see that it depends not only on contractions like (a A a B ) 

[cf. the second V diagram and the first T diagram], but also on (a A [3 B ac)o [ c f- the first 

V diagram, involving only scalar propagators]. Here again, we understand very simply why 

the 3-body interaction Lagrangian L B q of Eq. (|2.20cj ) depends only on these two types of 

A 

contractions, involved in G AB and (3 BC . 

The new parameters appearing at order G 3 can now be found by examining the diagrams 
connecting four bodies in Fig. [7|. While the H and F diagrams depend only on the previous 
contractions of oi AB and /3 AB , three new contractions occur in the first G, Z and X 
diagrams. Indeed, they involve respectively the contractions (/^a^a^a^o, (a a A (3^ b [3^a^)Q 
and (Robed® A a B a c a< b)Qi which are independent from (a A a B ) and ((y A p B ac)o m generic 
tensor-multi- scalar theories. It is convenient to introduce compact notations for these pa- 
rameters; we define 

£ a = ((3 l A a B a c aD)o , g * 

BCD [1 + (a A a B ) ] [1 + (a A a c )o) [1 + (a A a D ) ] ' 

Cabcd = (<xaPbPc<xd)q ^ ^ 

[1 + (<y A a B )o] [1 + («b«c)o] [1 + (otcot D )o\ 
Xabcd = (Rabcda A <y b B <y c c a D)o , (3.18c) 

where the choice of the carrier letters e, (, x is made by pictorial analogy with the cor- 
responding diagrams G, Z and X. The introduction of factors [1 + (aA«B)o] _1 is made to 
simplify some equations below, where we shall factorize the effective gravitational constant 



G AB , Eq. (pH). 



used to the notation, it is easier to see which new 2PN parameters can occur in the multi-scalar 
case rather than in the mono-scalar one. 
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The main conclusion of our diagrammatic analysis of tensor-multi-scalar theories is there- 
fore that only three new self-gravity- dependent parameters appear at the G 3 level. By con- 
trast, previous studies in the literature suggested the need for a much larger number of 
parameters in the general framework of Lorentz- invariant theories of gravity pl| , p2| . Let us 
note that if we had restricted ourselves to theories involving only one scalar field, only two 
self- gravity-dependent parameters e B Q D and (abcd would have showed up, the Riemann 
tensor of the scalar manifold being then identically zero 8 . We are going to see that the 
parameter Xabcd disappears anyway when considering the weak self- gravity limit. 



B. 2PK approximation 

The diagrammatic results of the previous subsection gave the structure of the non- 
linearity expansion, in powers of Gm/rc 2 , of the iV-body Lagrangian, keeping unexpanded 
all powers ofv/c and of the self-gravity s. The second post-Keplerian limit of tensor-multi- 
scalar theories can now be obtained by expanding the results above in powers of v 2 /c 2 (still 
keeping s unexpanded). First of all, the Minkowski proper time ta involved in ( |3.16| ) is 
obviously expanded as 



(Lta = dtA\ 1 — v^/c 2 = dtj 



2 c 2 



1 
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+ 



(3.19) 



The (v A /c 2 ) 3 term is necessary to write So[a] at the 2PK order, but the expansion ( |3.19| ) 
can be truncated to order (v A /c 2 ) 2 for the 2-body diagram "I", and to order (v\/c 2 ) for the 
3-body diagrams V and T. In the 2PK diagrams G, Z, F, H, X, it is enough to replace the 
proper time r A by the coordinate time t^. 

Similarly, one must expand in powers of v/c the unit 4- velocity u A = dx A /cd,TA 



Jc)/ a/1 — v^/c 2 , and in particular the contraction 2(f^ l/ u A u B ) 2 — 1 which appears for 



each graviton propagator connecting two bodies (see Appendix A). One gets easily 



2(u A u B f 



1 + 2 



+ 2 



(V4 - V B ) 2 K + V|) - (V A X V B ) 



+ 



(3.20) 



where x v B denotes the usual vector skew product. Finally, the Green function (|3.17|) 
can be expanded as 



cQ sym (x A -x B ) = [S(tj 

CO 

= E 

n=0 



|xa 



t B ~ tab/c) + S(t A - 
-* B {t B )\ 2n ~ l d 2n 



(2n)\c 



2n 



(dt B ) 



2n 



tB + r A B/c)]/2r AB 
8{t A - t B ) 



8 In the one scalar case, we could also express formally Qabcd m terms of 1PK parameters, using 
(aA/3,B/3c«D) = (oA/3Bac)(oB/3cQ ; D)/(oBOc)> but this expression is singular when as or ac — > 0. 
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S(t A -t B ) , \*a -x b (*b)| d 2 5(t A -t B ) 



+ 



r A s + 2c 2 



^1 



24c 4 



+ 



(3.21) 



The three terms in ( |3.21| ) are needed for the "I" diagram, while only the first two are needed 
for the V and T diagrams, and only the instantaneous Green function 5(t A — ^b)I t ab for the 
2PK diagrams G, Z, F, H, X. The time derivatives of 6 functions in ( |3.21| ) imply that the 
2PK Lagrangian depends not only on the positions x A and the velocities of the bodies, 
but also on their accelerations a^. (Higher derivatives can be eliminated, at 2PK order, by 
means of suitable integrations by parts.) As discussed in |T7 |, this is a consequence of our 
choice of gauge for the Einstein-Hilbert action S 8p i n 2- Indeed, all the accelerations can be 
eliminated by choosing for instance the "ADM" (Arnowitt-Deser-Misner) gauge [Hj, instead 
of the harmonic one we use in Appendix A below to simplify our calculations. (Actually, the 
simplest practical way of going from the harmonic gauge to a higher-derivatives-free gauge 
is to "wrongly" eliminate the higher derivatives by using the equations of motion in the 
Lagrangian [[T7|].) 

The 2PK Lagrangian describing the motion of N self-gravitating bodies is hence obtained 
straightforwardly from the action ( p. 14 ). We refer the reader to Appendix A below for its 
explicit derivation, and quote here only its general structure. It can be written as a sum of 
z-body interaction terms (1 < i < 4) 



L = L« + L® + L< 3 > + L< 4 > + 0(l/c 6 ) 



(3.22) 



where = ^2 A L^ A = — m A c 2 a/1 — v^/c 2 is the Lagrangian describing free bodies. 
Its 2PK expansion, generalizing ( p. 20a ) to order m A c 2 x (v A /c 2 ) 3 , is given by Q3.19 ). The 
2-body interaction Lagrangian L^ 2 \ given by the "I" diagrams of Fig. has the structure 



Gab t (2)g.r. 
~^~ l ab 



G AB m° A m B j AB ( (v A - v B ) 2 fi(r AB , v A ,v B ,a A , a B 



tab 



+ 



(3.23) 



where L AB G ' R ' is the expression obtained in general relativity (pure spin-2 interaction): 



Lab R = G*m° A m 



o ™0 

B 



1 f2(?AB, va, vb) f 3 (r AB ,v A ,v B ,a A ,a B ) 



tab 



c- 



c 1 



(3.24) 



and where the expression of the function f\ entering the term in Q3.23| ) is given in 

Eqs. ([AT7D-([AT8D. 

In other words, the 2-body interaction Lagrangian in tensor-multi-scalar theories 
presents two main differences with respect to the general relativistic result: (i) the bare 
gravitational constant G* is replaced by the effective one G A b, Eq. (|2.21a|) ; (ii) a correcting 
term proportional to l AB /c 2 must be added. 

The 3-body interaction Lagrangian, corresponding to the V and T diagrams of Fig. [7], 
has the structure 
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L (3) 



G AB G Ac m\m B wP c 



2 ^ 



1 + 2(3 



A 



TABTAC 



BC J , h PbcJa 



+E 



0{1ab,1aci1bc) 



A,B,C 



(3.25) 



where the functions / 4 and / 5 depend only on the positions and the velocities of the three 
bodies. The last term in Eq. fl3.25| ) denotes a sum of terms which are at least linear in the 
indicated ^ AB , an d which depend also on positions and velocities. 

Finally, the 4-body interaction Lagrangian, corresponding to the G, Z, F, H and X dia- 
grams of Fig. [7], has the structure 



L (4) 



\^ t (4)R.G. 1 \ ^ 
/ „ ABCD g / „ 



G AB G AC G AD m A m B m c m D 



A,B,C,D 



4 E 



A^(B,C,D) 

G AB GBcGcDm A m B m° c m D 

rjiBTBCTCDC 4 



r 'abt 'act 'adc 1 
(abcd 



'BCD 



24tt 



a^b^c^d 
1 x - Glm A m B mcm D 



c 1 



Xabcd / d x 



A,B,C,D 

OjlABilAC, 

C 4 
A,B,C,D 



*A, 



XcJ 



x — xa x — x# x — r x — xd 



+ E 



E 

A,B,C,D 



O ( AbC) ft ADi 



(3.26) 



where L AB ^?£ oc Glm A m B rriQm D is the general relativistic result, and where the 2PK 
parameters s B q D , (abcd and Xabcd have been defined in Eqs. ( pj.l8|) above. Note that the 
notation A 7^ B 7^ C 7^ D excludes only the equalities A = B, B = C or C = D. 



C. 2PN approximation 

The expression of the second post-Newtonian approximation can now be obtained by 
expanding the results of the previous subsection in powers of the compactnesses s ~ Gm/Rc 2 
of the bodies. To do this, we make use of results derived in |p^| . Following section 8 of this 
reference, we define 

c A = - 2 J- = 0(s) , (3.27a) 



c' A = 7^5 =0(s), (3.27b) 

a^-4^=0( S 2 ), (3.27c) 

cry 

b A = = 0(. 2 ). (3.27d) 
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Using the results (3.4) of |24| and neglecting the rotational kinetic energy terms with respect 
to the pressure (as is appropriate for solar system bodies), we can write these compactness 
parameters explicitly at 2PN order: 



c A = ^(U) A + - 4 (3^U-2f3U 2 ) +0(i) , (3.28a) 



c C \ <J / A \c 



o'a-^(U)a + o(- \ , (3.28b) 



c- \c 



= +o(- R ) . (3.28c) 



c 4 \a I A \& 



b A = - i (U 2 ) A + 0(-J . (3.28d) 



Here a = (T 00 + T u ) /c 2 denotes the mass density, p the pressure, U the potential sat- 
isfying AU = —4-nGa, and the angular brackets denote an average weighted by a : 
(/) = J afd 3 5t/ J ad 3 Sc. The tilde decorating the various quantities means that they are 
measured in physical units and expressed in terms of physically rescaled local coordinates 
x M = A((p(x A ))x^ adapted to describing the neighborhood of body A. The mass density 
a can be expressed in terms of the coordinate-conserved baryonic-rest-mass density, say 
Jm, as a — Jl + (Jth + 2p — JlU)/c 2 + 0(l/c 4 ), where h denotes the enthalpy. Note that 
Eq. ( p.28a| ) is more accurate than the usually employed 1PN expression for the compact- 
ness0: c A = (U) A /c 2 + 0(l/c 4 ) = -2E g A riw /m A c 2 + 0(l/c 4 ). A precise value of c' A is not 
needed for the following as it will only appear through the combination (4/5 — : y) 2 c' A . 

As shown in section 8 of [24]], the body-dependent parameters ( p. 21] ) can then be ex- 
panded as 

G AB /G = G AB /G = 1 - |(c A + c B ) + K + 4^ - |^ c A c B + P(2 + j)(a A + a B ) 

+ (| + C - 0) (b A + 6 B ) + 0(s 3 ) , (3.29a) 
Iab = 7 + V (l + | J (ca + Cb) + 0(s 2 ) , (3.29b) 
Pbc = P- {^ + \ + H-0 + Jh)c A - (^ + 0- f ) (cb + cc) 

-jc' A + 0(s 2 ) . (3.29c) 

Here we have set G = G*(l + afy for the Einstein-frame effective gravitational constant, 
and, as usual, r\ = 4(3 — 7. We have also introduced the notation 



9 Note also that it simplifies very much in the general relativistic case where the 1/c 4 correction 
vanishes. 
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These parameters are the weak-self-gravity limits of the parameters ( |3.18a|) -( [3.18b|) corre- 
sponding to the diagrams G and Z : 

e B £ D = e + 0(s) , (abcd = ( + 0(s) . (3.31) 

The parameter Xabcd, Eq. (|3.18cj ), vanishes in the weak- self-gravity limit because of the 
antisymmetry of the Riemann tensor R a bcd{f) — —Rbacdif) '■ 

Xabcd = (R abcd a a a b a c a d )o + 0{s) = + 0{s 2 ) . (3.32) 

Indeed, for the same reason, the first correction of order 0(s) is easily seen to vanish too. 

Our final conclusion is therefore that the 2PN limit of tensor-multi-scalar theories in- 
volves only the two new parameters e and (, Eqs. ( |3.30|) , besides the usual 1PN Eddington 



parameters (3 and 7. These parameters consistently enter into several 2PN effects. First, 
they parametrize two new independent contributions to the 4-body interaction Lagrangian: 

r W( c n- £ ST G 3 m A m B m° c m° D ( \^ G 3 m A m B m° c m° D 
b A^ { B,c,D) r Asr act ado 2 A ^ B ^ D r AB r BC r CD c 

Second, they parametrize the dependence upon self-gravity effects of the effective gravita- 
tional constant and of the 1PK parameters j AB and (3 BC . Discarding from Eqs. (|3.29| ) the 



2PN corrections proportional to the already experimentally constrained 1PN parameters (3 
and 7, we can rewrite them in the simpler form 

(7-,grav 7-igrav \ / zr-grav \ / pgrav 

tuaC m B c z J \mAC z J \m B c z 

| + C) {U2)a \ a {U2)b + 0(f3, 7) x 0(s 2 ) + 0(s 3 ) , (3.34a) 

Iab =7 + 0(^,7) x O(s) + 0(s 2 ) , (3.34b) 

Pbc =P + (e + C)^4 + C H4 + —2 + °(^T) x 0(a) + 0(s 2 ) . (3.34c) 
itiac 2 \m B c z m c c z J 



Eq. ( |3.34a| ) shows that e and ( are two independent 2PN generalizations of the well-known 



Nordtvedt 1PN modification of the gravitational coupling [second term of the right-hand 
side of ( p. 34a]) , r}(E A &v /rriAC 2 + E B ra,v /m B c 2 )}. In Eqs. Q3.34j ), we have simplified the writing 
by dropping all tildes in dimensionless ratios such as E A av /m^c 2 or (U 2 )a/c 4 . We shall 
continue doing so in the following each time this does not lead to any ambiguity. 

In fact the two roles (|3.33|) , (|3.34|) of e and ( are deeply connected. Indeed, the Lagrangian 



( |3.33| ) represents the interaction between an arbitrary number of non-self-gravitating mass 
points. From it we can formally reconstruct the Lagrangian representing the interaction 
between N weakly-self-gravitating bodies, A, B, . . . by considering each body as a collection 
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of mass points tua, ttib, ■ ■ ■ (held together by a slow orbital motion within the volume A). 
Let us denote by er(x) the average mass density within A, such that f A d 3 x.(r(x) = m° 

E 



A 



gives the total Einstein mass of body A. Then, if A and B label two point masses 
belonging to the same body A, one can rewrite the sum YIa^b Grn A m B /r AB as an integral 



AeA 



Ga (xx ) u (x 2 ) rf 3 X! <i 3 x 2 



Xi - x 2 



(3.35) 



where C/(x) = J rf 3 xG<r(x')/|x — x'| and where the angular brackets denote an aver- 
age weighted by o. Similarly, if A, B and C belong to the same body A, the sum 
J2a^b^c G 2 m A m° B m° c /rAB'rBC can be rewritten as an integral 



G 2 <T(xi)cr(x 2 )o-(x3)<i 3 Xi<i 3 X2<i 3 X3 



A 



o 



1 



, A |X X -x 2 | |x 2 -x 3 | 

Using such results, the interaction Lagrangian ( ^.33[ ) leads to 

G 3 m\m B m^m% ( 



(3.36) 



L(e,C) 



£ 

A^(B,C,V) 

+ E 

B^A^C 



G^m\m B m^m% e 

r AB r AC r AV C 4 6 



+ 



E 



A^B+C+V 



G 2 m A m B m^ 



E 



c 



3x £ - + C - 
6 2 



<t0 



r AB r BC r CT> c 



.4 



Gm A m B 



(U)a(U)b + 



3 e ( 
- x - + - 
2 6 2 



((f/ 2 )^+(f/ 2 ) B ) 



(3.37) 



where we have taken into account the different ways to choose infinitesimal elements in the 
same body to compute the correct multiplicities. [In pictorial language, if we denote the 
first two terms by (Ggcx>/3) and (Zabcv/2) respectively^ the following ones correspond 
to 3 x (e A g c /3) + {Zbaac/2), {Zbbac/2 + Z BACC /2), (Z AABB /2), § x + eJ A /3) + 

(Z AAA &/2 + Z A bbb/2), where a repeated index means infinitesimal volume elements inside 
the same body] 

The results ( |3.34| ) can be directly read off the Lagrangian ( |3.37| ). Indeed, the last sum in 

( p.37p leads to the e and ( renormalizations of the effective gravitational constant fl3.34a| ). 

while the 3-body sum J^b^a^c yields the renormalizations ( |3.34c| ) of the body-dependent 

A 

Eddington parameter (3 BC . 



IV. SECOND POST-NEWTONIAN ORDER METRIC 

A. iV-body physical metric 

As is well known, the iV-body Lagrangian contains enough information to derive the 
physical metric g^ at any spacetime point x x outside the bodies generating g^. Deriving 



10 The normalization of the S term is chosen for consistency with Appendix A. 
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this metric is important because it allows one to compute 2PN effects on, for instance, clock 
comparison, light deflection or time-delay experiments. To compute g M 2,(x A ), one introduces 
a test particle of negligible mass^] rfiQ located at x x = (x°,x l ), with an arbitrary 4- velocity 
= (c,v % ). We can now write the iV+ 1-body Lagrangian describing the N massive bodies 
and this test particle, and identify its mo-dependent part with the individual Lagrangian of 
the latter, namely 



L tcst - particlc (x\t/) = -fhocJ-g^ix^v^ 

= -m c 2 [-goo(x x ) - 2g oi {x x y/c - ^(z^V'/c 2 ] 1 . (4.1) 

This identification therefore allows us to compute 'g^ v {x x ) up to order 1/c 6 included for g Q0 
[i.e. up to 0{Gm/rc 2 ) 3 included], 1/c 5 for g 0i , and 1/c 4 for g^. 

If we are concerned only by the dependence upon the new parameters e and ( appearing 
at the 2PN order, Eq. ( |3.37| ), we can consider a particle at rest, since the e and ( terms are 
not velocity dependent. We therefore conclude immediately that g 0i and g^ do not involve 
the parameters e and ( at the 2PN approximation (i.e., at order 1/c 5 for goi and 1/c 4 for 
gij), and that goo = — 1 + 8 goo can be deduced from the Lagrangian 

^est-particlc^A) = _~ qQ 2 ^ _ ^ ^ = + ^fh C 2 ^ 00 + O(6g 00 ) 2 . (4.2) 

Comparing (fO) with (^.37[) yields 



x ~ i s\ £ G 3 m A m B m c , \- G 3 m A m B m c , _ sr^ G 2 m A m B 

5 2 PN^00(^, C = o 6" + £ 1^ 6 + 2£ 1^ 6" \ U A 

3 r oAroBTo C c & A ^ c) r 0A r AB r Ac c* ^ B r 0A r ABC s 

+2( \^ G 3 m A m B m c + G 3 m A m B m c 
A he r oAro B r BC c 6 r 0A r AB r BC c^ 

+c E ^%mA + (v)b) + 2c E ^^((^ + (^b) 

+ (s + 2C)J2^(U 2 ) A , (4.3) 

where r 0j 4 = |x — x^|. [All the masses entering this equation are evaluated at ip = <po, 
although we have dropped their index o for easier reading.] When a continuous description 
of the source bodies is used, this takes the simpler form 

Ww(£ , = _L,3 (X) + ^ + | j„ j^^fl 

,, 4) 



11 The index o of this test mass should not be confused with the time component x°. It has been 
chosen so that the body-dependent quantities a a A , ft a A h , . . . of Eqs. ( 2.17 ) reduce to their background 
a a ((po), Po b = /3 ab (ipo), ... for the test mass mo- 



values «q 
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Using a self-explanatory notation, the four terms of ( f4.4|) correspond to the diagrams E AAA , 
E AA0 , 7-aaao and Z AA0A respectively, the index a meaning here any infinitesimal volume of 
matter, and the index o referring to the spacetime point x x where the metric is computed. 
It should be noted that formulae (|4.3|) and ( |4.4|) are valid in the g*-harmonic gauge used 
to derive the iV-body Lagrangian in subsection III-A. However, it is easy to see that the 
corrections proportional to e and ( are the same in all usual coordinate systems (ADM- 
gauge, g-harmonic, . . . ), since these parameters appear for the first time at order 1/c 6 in the 
time-time component of the metric^. In particular, the fact that (fti an d gij do not involve 
e and ( at the 2PN order is valid in any of these coordinate systems. 

Let us mention for completeness that the presence of e and ( modifies the total energy 
of a gravitating system by the amount 

feN%,0 = ~ f - ± |/^? ggg|g|gg) , (4.5) 

where the tilde decorating the various quantities means as before that they are measured in 
the physical units (|2.6|) . 



B. Exact one-body metric 



In the present subsection, we verify and complement the above results by using the 
exact static and spherically symmetric solution of the field equations (|2.8a|) - ([2.8b| ) derived 
in section 2 of f23| . In the coordinate system introduced by Just P^ , |SD[| , 

dsl = -e v c 2 dt 2 + e~ v [dr 2 + e\d9 2 + sin 2 9 dcf) 2 )] , (4.6) 

we found the solution 

e A = r 2 - ar , (4.7a) 

e» = (1 - a/r) b ' a , (4.7b) 



where 

b - = ^ / v^ 3 x(-n + 2i) (4.8) 



is the Einstein-frame mass parameter, and a is a constant of integration. Introducing the 
parameter p = (1/a) ln(l — a/r), the equations controlling the radial variation of the scalar 
fields can be expressed by saying that <p a (p) follows a geodesic of the metric jabif ), i-e., 
5 J , y a } ) ((p c (p))(d(f a /dp)(dif b /dp)dp = 0. In particular, p being an affine parameter, the norm 
of d(p a /dp is constant and given by 



Of course, it is always possible to introduce by hand a spurious dependence on these parameters, 
by redefining for instance the spatial variables as x H = x l +0{e, C)/c 4 ) whereas any usual coordinate 
transformation involves e and £ at order 1/c 6 . 
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±lab^)%-%- = a 2 -b 2 = const . (4.9) 
dp dp 

The actual geodesic traced out by ip a (p) in the internal scalar space is uniquely determined 
by the initial values at p = (i.e., at spatial infinity): f a (p = 0) = ip a (r = oo) = y?g and 
(d(p a /dp)(p = 0) = k a , say. Assuming and k a given, we can choose field variables in the 
scalar space so that the metric 7ab(v 9 ) reduces to the constant 6 ab all along the line v 5 ( r )re[o,oo] 
(Fermi coordinate system in the 93-space). This allows us to write the solution of (|4.9| ) in 
the simple form 

(p a = <p a +pk a = tp a + — In fl - -) . (4.10) 



Using then the results of |2~4jl , we can relate the n integration constants k a to the matter 
distribution: 

k a = % / ^ 3 xa»(-T° - T;j . (4.11) 



Note that ( |4.11|) is valid only if the (^"-coordinates are "Fermi" all along the solution (p a (r). 



When using generic scalar variables, the quantity a a ((f) appearing on the right-hand side 
of (|4.11|) must be replaced by the parallel-transported value of the vector a a (<p) up to the 



point (fQ. To simplify, the reader can think that we work from the beginning with a flat 
cr-model metric 7 a fe(v 9 ) = <W The only 2PN contribution that we would forget in that 
case would derive from the first X diagram of Fig. [j], and would involve the contraction 
Rabcd(H ) ) a \ ab A aC A a o which vanishes identically because of the antisymmetry of the Riemann 
curvature tensor R abcd (ip) = -R b acd{v)- 

By comparing the behavior at infinity of the solution ( |4 . 1 0| ) 

^• = ^-£ + 0(1) (4.12) 

with the one obtained for an isolated pointlike body A (see Eq. (|2.18|) ) 



we can deduce that k a = G^m^a^/c 2 so that, using ( [4.8| ) and ( [4.11| ), the actual expression 
of the parameter a a A for extended bodies reads: 

„ _ 2k a _ /V^d 3 xa a (y>)(-T°-7^) 

b ivrJM-n + Ti,) ■ l4 - i4j 

Note the change of sign of the pressure term T\ { between the numerator and the de- 
nominator, which makes a a A differ from at order 1/c 2 even if the function a a (ip) is 
a constant, as in the Jordan-Fierz-Brans-Dicke theory. In the latter case, one finds 
a a = a a x (1 — 2 J ^fg^d^xTlJiriAC 2 ). When comparing this result with = a° a + 
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(dm A /d\nG)(d\nG /d(f a ) = a a x (1 — c A ), we get the following expression for the com- 
pactness in Jordan-Fierz-Brans-Dicke theory: 



ca = 2 J y^d^TlJrriAC 2 . 



(4.15) 



This result is an (exact) generalization of the Newtonian virial theorem J c/ 3 xcr[/ = 
6 J <i 3 xp + 0(1/ c 2 ) which is valid in the Jordan-Fierz-Brans-Dicke theory and therefore 
also in general relativity (in the limit a a — > 0). 

The exact solution (f4.7|) -( 4TT0[) can now be expanded a la Eddington in powers of 1/r. 
Although the coordinate system ( f4.6|) is as good as any other to verify the 2PN terms 
involving the new parameters e and (, it is convenient to express our results in isotropic 
coordinates, to make better contact with the general relativistic result. Let us then define 
a new radial coordinate p such that p(l + a/4p) 2 = r. The Einstein line element (|4.6) now 



reads |$9| 



ds* 



1 



4p) 



\ 2b/a 



1 + 



-2b /a 



c 2 dt 2 



a 

Jp 



2-26/a 



+ I 1 

-c 2 dt 2 
+ [dp 2 + p 2 (d6 2 + sin 



4p 

\ 2+26/a 

l + Tp) [dp 2 + P 2 (d6 2 + sin 2 



! b b 2 b 3 
p 2p 2 6p 3 



l + a 2 A 



b 

1 + - 

P 



2p 2 



] )] 



l + a 



(4.16) 



where we have used Eqs. ( |4.9|) and (|4.14j ) to replace the constant a 2 in terms of b 2 and 



a 2 A = (a^7 a ba^) ¥ ,= ¥ , . The metric which determines the dynamics of test particles in the 
vicinity of A is not the Einstein metric g* but the physical one g^ u = A 2 ((p)g^ u . We must 
therefore expand also A 2 ((p) in powers of 1/p, by using the exact solution fl4.10| ) of the scalar 
fields. We get 



A 2 (<p) 
A%o) 



1 - (a A a )- 
P 



{a A a ) 2 + 



(a A f3 a A ) 



2p 2 



(Paaa) + G(a A a )(a A p a A ) + 4(d A d ) 3 + 



(audio) (1 + a 



b 3 
24P 1 



+ 



(4.17) 



where we have set (d^do) = (d^7 ab d% =¥ , , (a A f3 a A ) = (a A f3 ab a b A ) v=Vo , (!3' abc a A a A a A ) v = Vo , 
and where the a-model tensors a a , [3 a b and f3' abc have been defined in Eqs. (|2.12|) above. 
Finally, when using the physical coordinates p = A(Lp )p and t = A(<p )t (see Eq. (p.6[) ), we 
find that the physical line element ds 2 = A 2 ((p)ds 2 reads 



d's 2 = —c z dt 



2j72 



, Pa 

pc z 



pc z 



B) 



27 



+ [dp 2 + p 2 (d6 2 + sin 2 



x 



20(1 + 1 M ) 



Pa 
pc 2 



;i + r) + o 



where 

B = 9^ + 3 Paa + (4 + 7ao) ^ " (1+W2) "l8~ 

r = /aa + (28 + i5 lA0 )— + (1 + W2) -y , 

and where we have set 

AU = ha/Mvo) = G A0 m A = G*A 2 (ip ) [1 + (a A a )] m A , 



7ao = -2 



1 + (a A a ) 



Iaa = -2 



a. 



l + a 2 ' 



O 1 2 

Paa = -^{ot A (3 a A )/ [1 + (a A oio)] 
o 

AAA 



(/3lb c «A«A«A)/ [1 + («A«o)] 1 



(4.18) 

(4.19a) 
(4.19b) 

(4.20a) 
(4.20b) 

(4.20c) 
(4.20d) 



As a check of the consistency of our method, it is useful to compare the 2PK Lagrangian 
of subsection III B with the test-particle Lagrangian (|4.1|) written for the metric QL18Q . Let 
us first rewrite this metric in g*-harmonic coordinates. We find 



-^oo = l-2^4 + 2 

rc \rc / 



(re 2 ) 



3B 



(l + Paa 
7ao( 4 + 7ao) 



+ 



7aa 



'(2 + 7aa) 



9ij 



x 



l + 2^( 1+ ^) + (M) 2 



4(2 + l AA/ 
o 



+ 
7 



1 

=7 \2 



i + 2p AA + 3j A0 + ~tf A0 y + 



(4.21a) 

7aa(2 + 7ao) 2 



- 2 vrcV 2(2 + 7 AA ) V C V 



4(2 + 7aa) 

(4.21b) 



where r = a/ <5ijX l xi denotes the g*-harmonic radius, related to the isotropic radius p by 

Ha 



P 



1 + 



2pc 2 



1+7aa/2 



(4.22) 



[This g*-harmonic radius r should not be confused with the Just radius used in Eqs. Q4.6Q — 
( |4. above, also denoted by r.] The test-particle Lagrangian deduced from the metric 
( 4.21] ) agrees with the Lagrangian (|3.22| ) — ( |3.26|) , written in the particular case of two bodies 
m , m A . It is notably easy to check the I^^a^ contribution in (|3.26|) , which comes from 
the |B term of ([£2IiiD . [The Caoao and 

Xaoao contributions in (|3.26|) cannot be checked, 
because they involve the square of the test mass m .} To ease the reading, we have used a 
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slightly inconsistent notation in expressing the metric coefficient g^ (corresponding to the 
coordinate system x) in terms of the original (Einstein-frame) coordinates x. We use here, 
as we did above, the simplifying fact that dimensionless ratios such as /i/rc 2 are numerically 
equal to their physical-units counterparts Jl/rc 2 . 

The 2PN (weak-self-gravity) limit of the one-body metric ( |4.18| ) [or ( [4 . 2 1| ) ] can be ob- 
tained by expanding the body-dependent parameters in powers of the compactness of body 
A, as in Eqs. ( |3.29| ) above: 



@ AA 



B 

7ao 



Gm, 



l-|c A + /?(2 + 7)a A +(! + C 



i/r & 



+ 0(s s 



(3-(( + 2f3-r ] f3)c A + 0(s 2 ) 



& H 2 A (1 + (3 AA ) = (Gm A ) 2 (1 + $ - (C + 6/3 - 7)ca + 0(s 2 )) 



-e + -(3 + — +0(s) , 
9 3 M 18 K J 

7 + 77(l + |) C A + 0(S 2 ) 



IIa(1+Jao) =Gm A (1 + 7 + |c j4 + 0(s 2 
15 + 87_ 



P 



6 



7 + 0(s) . 



(4.23a) 

(4.23b) 
(4.23c) 

(4.23d) 
(4.23e) 



We thus recover that the spatial metric g^ does not depend on the parameters e and ( 
at the 2PN order, as shown by Eqs. ( ^4.23d| ) and ( |4.23ej ). On the other hand, the e and ( 
contributions to g 00 can easily be deduced from the above results, 



52PN0oo(e,C) 



(e + 2C) Gm A 



-C(/3,7) 



(U 2 )a + 



2C {Gm, 



P 



Gm, 



(4.24) 



consistently with Eqs. 



V. 2PN EXPERIMENTS 

A. Constraints from 1PN experiments 

Solar-system experiments impose tight bounds on the Eddington parameters 7 and j3 
| 28| , |29[| . Using Eqs. (|2.15b| ), they can be written as 



|7| < 2 x 1(T 3 -IO- 3 < -7—^ = (a a -? ab a b ) < 10" 3 , (5.1a) 

2 + 7 

\(3\ < 6 x 1(T 4 -1.2 x 1(T 3 < - — = (a a (3 ab a\ < 1.2 x 1(T 3 . (5.1b) 

(2 + 7) 



Let us first discuss the lessons one can draw from the bounds ( |5.1| ) if one interprets them 
within the most natural theoretical framework. Since we assume that the (a-model) metric 
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7afe is positive definite, so that the theory does not incorporate any negative-energy excita- 
tion, Eq. ( 5. lap constrains the magnitude of all the components of ag. 1PN experiments 



thereby constrain the linear interaction of the scalar fields to matter to be small. On the 
other hand, our diagrammatic analysis of section III shows that any observable deviation 
from general relativity involves at least two factors a^, to fill the end blobs connected to 
scalar propagators, in diagrams such as Fig.|7| or any higher order ones. Moreover, theoretical 
"naturalness" leads us to expect the coupling function A(<p) entering Eq. ( ^.3[ ) not to involve 
any large dimensionless number, so that the successive derivatives of In A{tp) [such as /3° 6 , 
P'abci ■ ■ ■) are a V™ r i expected to be of order unity. The conclusion is that we anticipate 
the new parameters entering any post-Newtonian order to be at most of the same order of 
magnitude as 7. In particular, the 2PN parameters ( |3.30| ), which involve respectively 3 and 
2 factors «q, are expected to be of order |e| ~ | ct 1 3 ~ 3 x 1CT 5 and \(\ ~ «q < 10~ 3 . The 
corresponding 2PN deviations from general relativity are thus generically constrained at the 
level ^5 10~ 3 x (Gm/Rc 2 ) 2 ~ 1CT 15 , too small to be detectable in the solar system, even 
in future high-precision experiments. From this point of view, our conclusion is therefore 
that solar-system tests cannot probe the 2PN structure of gravity. However, they can give 
a clean access to the Eddington parameter 7, which is of greatest significance as it measures 
the basic coupling strength between matter and the scalar fields, and constrains all the other 
PN parameters. Within this viewpoint, the only meaningful testing ground for 2PN and 
higher-order effects are binary-pulsar experiments. Indeed, since the self-energy Gm/Rc 2 
of a neutron star is typically of order 0.2, the 2PN effects ^ 10~ 3 x (Gm/Rc 2 ) 2 can yield 
significant deviations on timing data. Moreover, in such strong-field conditions, the sum of 
the series of all higher post-Newtonian orders may be large enough to compensate the global 
«o factor of all the parameters. We have indeed shown in [4(| that non-perturbative strong- 



field effects can compensate even a vanishingly small «q ~ — ~j. The overall conclusion is 
that two different regimes of gravity must be distinguished: (i) the weak-field conditions 
of the solar system can give a clean access to the fundamental parameter 7; (ii) higher 
post-Newtonian effects can be tested in the strong-field regime of binary pulsars. 

These conclusions, derived from a theory-motivated point of view, happen to be valid 
if we adopt a more phenomenological point of view, i.e., if the coupling function A((p) is 
now supposed to have a priori any possible shape. It can thus involve in principle large 
dimensionless numbers, and the 2PN parameters ( |3.30| ) can be of order unity in spite of the 
bounds ( p. la]) on 7. This is particularly true of the parameter (, which can be of order unity 
even if A(ip) does not involve numbers larger thanQ ~ 30. For instance, in a model with 
two scalar fields and a flat cr-model metric ■y a b = S a b, the values 

al^W- 3 , a 2 = , fti «1 , A 2 ^30 , (5.2) 

give (a 2 )o = ol\ rs 10~ 3 and (a/3a)o = ai/?ii«i ~ 10~ 3 , consistently with (|5.1[) , but ( ~ 
(a/3/3a)o = ai/^ai + a\(5\ 2 a\ ~ 1. By contrast, getting e ~ 1 seems less natural, as some 
components of /3'® bc would have to be as large as 3 x 10 4 to compensate the «o«o a o f ac t° r - 



13 Let us note in this respect that in the string-derived model of Ref. |26|, the equivalent of Pu is 
expected to be of order 40. 
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More generally, it is easy to see that ( is a priori unconstrained by the solar-system 



bounds (|5.1|) . Indeed, from the positive-definiteness of the metric 7^, we can use the Cauchy- 



Schwarz inequality to obtain a lower bouncQ for £: 



(a(3(3a) > [(a(3a) Y K & (> 8(3 . (5.3) 



^2 



On the other hand, ( has a priori no upper bound, and we can have ( > /I7I if the 
(y9-space) vector (/3 a b a b ) = (/3a)o is almost orthogonal to «q, as in the example ( p.2| ) above. 



In Ref. we adopted an even more "phenomenological" point of view, by assuming 
that the a-model metric 7^ could have a priori any signature (i.e., that the theory could 
involve negative-energy scalar fields). In that case, the bounds (|5.1a|) on 7 do not constrain 



the magnitude of Oq, but only its direction: It should be almost a null vector in the yj-space. 
Then, it is easy to see that e is not constrained at all by the experimental limits ( |5.1| ), and 
that ( is not constrained either if the theory involves at least three scalar fields. However, 
from a theoretical point of view, it seems more plausible that the gravitational interaction 
involves only a small number of massless scalar fields. If it involves only two of them (or 
less), the Cayley-Hamilton theorem can be written as 

(af3 2 a) = (Tr/3)(a(3a) - (det/3)(a 2 ) , (5.4) 

where the trace and determinant are taken for the 2x2 matrix (f3 a b ). This relation shows 
that even in the case of a hyperbolic metric 7^, the parameter ( is constrained to vanish 
with f3 and 7. 

The conclusion of the above discussion is therefore that the 1PN experimental limits 
( p.l[ ) most plausibly constrain the 2PN parameters e and ( to be very small, but that 
there is no theoretical impossibility that the latter be of order unity. In the following, 
we will adopt a phenomenological point of view, and consider that these parameters are a 
priori unconstrained to discuss the possible experiments which may measure them. Let us 
underline, however, that this "phenomenological" attitude is still in the framework of our 
field-theory approach, and should not be confused with previous studies in the literature 
l2| which aimed at describing any conceivable deviation from general relativity. 



B. Solar-system experiments 

1. Light deflection and time delay 

As discussed in the previous subsection, the most important parameter to determine 
experimentally is 7, as it measures the basic coupling strength and constrains a priori 



all the possible deviations from general relativity. It has been shown in ||27|j26|1 that the 



cosmological evolution of tensor-multi-scalar theories generically drives 7 towards 0, and 



The fact that we are able to derive relations between 1PN and 2PN parameters is a further 
illustration of the power of our field-theory approach. 
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that its present value is expected to be ^ 1CT 7 . The 1PN deviations proportional to 7, 
entering light-deflection and time-delay experiments, are thus expected to be comparable 
to 2PN effects, if we adopt the phenomenological point of view that the 2PN parameters e 
and ( are a priori unconstrained (and that j3 may be much larger than 7). It is therefore 
important to determine whether these effects can be distinguished from each other. We shall 
see here that light-deflection and time-delay experiments are ideally suited to an accurate 
determination of 7 in the solar system. 

First of all, it is immediate to see that these experiments do not depend on e and ( 
at the 2PN level. Indeed, our results of section III and IV show that these parameters 
appear for the first time at order 0(l/c 6 ), in the time-time component of the metric goo • 
Therefore, they do not enter the 2PN geodesic equation satisfied by light, which involves 
only the 0(l/c 2 ) and 0(l/c 4 ) orders of the metric. Consequently, any experiment probing 
the propagation of light is independent from e and ( at the 2PN level. 

We compute in Appendix C the 2PN light deflection by a massive body and find 



A<f> 



2/^(2 + 7^) 
Poc 2 



/M2 + 7) 



PoC z 



+ 



7T 



Poc 2 



2 r 



15 + y 7 + 47 2 



O 



(5.5) 



Here po denotes the minimal distance between the light ray and the center of body A 
measured in isotropic coordinates. The body-dependent parameter j A0 entering the first 
term of ( |5.5[ ) can be expanded in powers of the self-gravity of body A, as in Eq. ( |4.28d| ) 
above: 



Iao = 7 - (4/3 



7)(2 + 7)^r V MAC 2 + 0(l/c 4 ) . 



(5.6) 



Note that /3 appears in the self-energy corrections to the 1PN effects, whereas it is absent 
from the 0(Gm/rc 2 ) 2 contributions. 

Contrary to the hopes of Refs. [|l2|,[r|,|l^,[n| , a conclusion of our approach is that improved 
light-deflection (or time-delay) experiments do not give access to any theoretically significant 
2PN parameter. In fact, the formal 2PN generalization of the Eddington 7 parameter 
introduced in Refs. P^JT3[ ] (under the different notations e and 5; see also [15] where it is 
denoted A) is equal to 1 + T, where T is the function of 1PN parameters given in Eq. ( |4.23c|) 
above: 



^Epstein— Shapiro ^Fischbach— Freeman 



A 



. 4- 15 + 87_ 

Richter-Matzner — 1 + ~P H ~ 7 ■ (5.7) 

3 



Our conclusion concerning the impossibility to probe significant 2PN deviations from 
general relativity in light-deflection and time-delay experiments should not be viewed only 
negatively. The positive aspect is that these experiments can give a clean access to the 
fundamental parameter 7, i.e., that no 2PN deviation from general relativity can complicate 
their interpretation in the case where the data are too scarce to allow a clean separation 
of 1/r and 1/r 2 effects. Indeed, at the level 7 ^$ 1CT 6 that these experiments aim for, the 
2PN effects proportional to this parameter in Eq. (|5.5| ) are totally irrelevant^. On the other 



15 



At the qualitative level, it remains that checking the coefficient 15ir/A appearing in the second 
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hand, one remarks that the 7 parameter which will be measured by these experiments is the 
body-dependent quantity 7 A0 rather than its weak-field limit 7. However, the self-gravity 
renormalization of 7 displayed in Eq. ( |5.6|) is already strongly constrained by the LLR (Lunar 
Laser Ranging) bounds on 77 = 4/3 -7 (\rj\ < 2 x lCT 3 fl29fl ). Using |£| rav |/m c 2 ~2x lCT 6 
for the Sun, one gets 

(4^-7)(2 + 7)l^ raV |MoC 2 < 10- 8 • (5.8) 

Therefore, an experimental determination of 7 0O at the 10~ 6 or 10~ 7 level, which is expected 
to be reachable in missions presently considered by the European Space Agency such as 
GAIA (Global Astrometric Interferometer for Astrophysics) or SORT (Solar Orbit Relativity 
Test), will indeed give a clean measurement of 7. 

Finally, let us remark that the value of 7 we are talking about in this work is the 
one corresponding to the value of the scalar-field background <£>q around the solar system. 
As remarked in [p2|, it differs from, say, the cosmological average of 7 (which is the one 



discussed in Refs. |2^J2^j) by the effect of the spatial fluctuation AU of the gravitational 
potential. However, we show in Appendix B that the corresponding change in 7 is given by 
4/3(2 + 7) AU/c 2 and is therefore expected to represent only a very small fractional change 
of 7 of ordeiQ \/3 a b \AU/c 2 , where \/3 a b \ is the norm of the matrix f3 a b (say the modulus of its 
largest eigenvalue). 



2. Other possible 2PN experiments 

The previous subsection showed that experiments on the propagation of light in the 
solar system cannot give access to any post-Newtonian parameter but 7. The question that 
we will address now is whether the 2PN parameters e and ( can be measured in the solar 
system, or at least if it is possible to constrain their magnitudes at an interesting level. 

One of the most famous tests of post-Newtonian gravity is the perihelion shift of Mercury, 
and more generally tests obtained through a global fit to the orbital motion of the planets. 
Before any calculation, it is clear that we cannot hope to measure the Mercury perihelion 
advance by traditional means to an accuracy of 2PN significance. Indeed, the 2PN effect is 
smaller than the 1PN prediction by a factor ~ (Gm Q /rc 2 ), and corresponds to an advance 
of ~ 10~ 6 arcsecond per century. However, it may be, one day, possible to reach this level 
by using an artificial satellite orbiting around Mercury and tracked with very high accuracy 
from Earth. To increase the parameter Gm Q /rc 2 , one could also observe the perihelion shift 
of a drag-free satellite in close elliptical orbit around the Sun, but the construction of such 
a satellite seems unrealistic with present technology. Anyway, several reasons show that 
such experiments could not constrain e and £ at an interesting level. Indeed, we compute 



term of ( |5.5| ) will be a new confirmation of the nonlinear structure of general relativity, even if it 
does not constrain any plausible theoretical alternatives. 



16 Even if \P a b \ ~ 30, the fact that AU cosmo /c 2 ~ 10 5 indicates a negligible fractional change. 
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in Appendix C below the perihelion shift per orbit for a test mass mo <C m®, and we get in 
isotropic coordinates 



A0 6nGm& 



a(l — e 2 )c 2 



27 - + Cc© Gm /7 £ , 2 , 
1 + 3 + ^ 2 + 6 +0(e ' 



c 4 \ c u 



(5.9) 



where a is the coordinate semi-major axis of the orbit, e its coordinate eccentricity (in 
isotropic coordinates), and c ~ 4 x 1CT 6 is the compactness of the Sun. Although ( enters 
this expression, it cannot be distinguished from a small contribution of 27 — (3 at the 1PN 
level. In fact, it comes from the expansion ( 4.23b ) of the body-dependent parameter /3 00 = 



(3 — (c & + 0((3, 7)/c 2 + 0(l/c 4 ), and is thus a mere renormalization of (3. On the contrary, the 
contribution proportional to e can in principle be distinguished from the 1PN effects, since 
it involves a different power of (Gm/a). Numerically, we find that Mercury's perihelion 
is deviated from the general relativistic prediction by ~ 0.5e mm per year, too small to 
be of observational significance. In the totally unrealistic situation of a drag-free satellite 
grazing the surface of the Sun, this 2PN deviation would be of order rs 30e m.yr -1 , which 
is a priori much easier to detect. However, to distinguish the 0(Gm/ac 2 ) and 0(Gm/ac 2 ) 2 
contributions, it would be necessary to compare several satellites at different distances from 
the Sun, or alternatively to look at periodic effects on a given orbit, which may be much 
more difficult to observe than secular effects. Moreover, one will have the difficulty of 
separating the 2PN contribution from the Newtonian contribution due to the quadrupole 
moment of the Sun, which as the same dependence^ on a (not to mention the huge thermal 
and electromagnetic effects of the Sun on a close satellite). In conclusion, perihelion-shift 
experiments in the solar system can in principle give access to the 2PN parameter e, but 
present technology does not give the hope of measuring it to any significant level. 

We have seen in Eqs. ( |3.29a| ) and ( |3.34a|) above that the effective gravitational constant 
Gab depends on the self-energies of the bodies A and B. In particular, the Earth (©) and 
the Moon (C) do not fall with the same acceleration towards the Sun, since G ® 7^ G Q c- 
This violation of the strong equivalence principle implies a polarization of the Moon's orbit 
(Nordtvedt effect), that can be tested in the Lunar Laser Ranging data. The deviations 
from general relativity are proportional to the ratio 

G&(B g G&C = ~|( c ® ~ c c) + (C + 0(P, 7 y)) c (c© - cc) + 0(a©,6®,a c ,&c) , (5.10) 

where cc < c$ <C c are the compactnesses of the three bodies, and where we have neglected 
<% ~ fr® ~ c | ( as wen as a c an d be) with respect to c c . The first term, involving the 
parameter 77 = 4/3 — 7, is the standard 1PN deviation which has been constrained at the level 
\tj\ < 2 x 10~ 3 by LLR data. The dominant 2PN contribution, C c ©( c ® — c c)> is equivalent 
to a renormalization of 77 into (77 — 2(c Q ) ~ (77 — 1CT 5 C). From a theory-based viewpoint, 
such a renormalization is of no consequence: as rj and ( are both proportional to the basic 



17 The peculiar anisotropy of multipolar effects will, however, help in this respect. 



34 



coupling strength 7, we can always neglect the fractionally small correction 10~ 5 £ to 77 (in 
the same way, as explained above, that one does not have to worry about "cosmic variance" 
effects). However, from a phenomenological viewpoint, ( is an independent quantity which 
would complicate the interpretation of a high-precision LLR experiment reaching the 77 ~ 
10~ 5 level. [This level would correspond to measuring the Earth-Moon distance with 0.1 
millimeter accuracy] In other words, the phenomenological point of view obliges us to look 
for other independent experiments allowing one to separate the effects of rj and (. 

The polarization of Mercury's orbit around the Sun due to the presence of Jupiter (i.e., 
the Nordtvedt effect for Mercury) can give access to a different combination of rj and the 
2PN parameters e, (. The corresponding deviation from general relativity is proportional 
to 

Gjo — Gjm V (£ 



G 



■|c + (| + () b G + C x 0(c qCj ) + 0(j3rf)0(c M , c|) , (5.11) 

where cm <C cj <C c & are respectively the compactnesses of Mercury, Jupiter and the Sun, 
and where we have neglected c Q cj with respect to 6 Q ~ cj|. The dominant 2PN contribution 
plays again the role of a renormalization of 77, but this time into [77 — (e + 2C)6 Q /c Q ] ~ 
[rj — 10~ 5 (£:/2 + £)]. Note that e enters now this expression, whereas it was absent from the 
corresponding lunar result. This complicates the problem of separating the contributions 
of 77 and (. It would be necessary to dispose of a third experiment giving access to e, for 
instance by using the perihelion shift (|5.9|) . However, we have seen that such a measure is 
not likely to be performed with current (or foreseeable) technology. 

Another attempt to determine e could be to compare ultra-stable clocks: one located on 
Earth and another one somewhere close to the Sun. Indeed, the Einstein effect gives access 
t° a/— S'oO) where —(Too is given by Eq. ([4. 2 la]) in g*-harmonic coordinates and by the first 



bracket of Eq. ( 4.18[ ) in isotropic coordinates. The general relativistic prediction for the rates 
of clocks is thus multiplied by a factor [1 + (Gm Q /rc 2 ) 2 /3 QQ — ^(Gm Q /rc 2 ) 3 e + 0(0, 7)/c 6 ]. 
This gives (relative to a clock on Earth) effects of order « 1.6 x 10 _18 e for a clock located at 
the surface of the Sun, and ~ 3 x 10~ 2A e for a clock on Mercury. With foreseeable technology 
(10~ 18 stability), one might barely be able to constrain e at the 0(1) level. As for the 
perihelion shift ( |5.9|) , it should be noted that the parameter ( contained in /3 00 ~ (3 — (c & 



is a mere renormalization of j3 and cannot be accessed independently. 

A possible way to access this ( contribution in j3 AA would be to measure both j3 Q& ~ 

(3 — (c & and /3^ e ~ (3 — £c© ~ (3 at a level < c ~ 4 x 10~ 6 . This would necessitate the 
tracking of Mercury with few cm accuracy, and the observation of a drag-free artificial Earth 
satellite at the 10 -2 cm level ||41|| . This second condition is two orders of magnitude smaller 
than what can be done presently. 

The conclusion of the above discussion is that the 2PN parameters e and ( are extremely 
difficult to measure in the solar system. Within a phenomenological approach, the only role 
of these parameters is negative: they complicate the interpretation of high-precision 1PN 
experiments. Hence the solar system appears not to be an appropriate testing ground for 
probing the 2PN structure of gravity. On the contrary, it is perfectly suited for measuring 
the fundamental parameter 7, as underlined in subsection V-B-l. 
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C. Binary-pulsar experiments 



Let us now consider binary-pulsar experiments, which will turn out to be much better 
testing grounds for the 2PN structure of gravity. Since the self-energy of a neutron star 
is typically of order 0.2 (as compared to ~ 10 -6 for the Sun), the 2PN and higher-order 
effects play an important role in the behavior of binary pulsars. In the present subsection, 
we will show that the combined analysis of several binary-pulsar data has the capability of 
constraining both e and ( at an interesting levelQ. Our aim is not to perform a full statistical 
analysis, but rather to illustrate the different types of constraints that can be obtained. We 
will consider four binary pulsars, for which different observable quantities can be measured. 

The Hulse-Taylor binary pulsar PSR 1913+16 has been continuously observed since 
its discovery in 1974. Besides the Keplerian orbital parameters P (orbital period) and e 
(eccentricity), three "post-Keplerian" |3l| , |32]| observables have been measured with great 
accuracy: the periastron advance u, the secular change of the orbital period P, and the 
time-dilation parameter 7timm g which describes both the second-order Doppler effect due to 
the velocity of the pulsar and the redshift due to the presence of its companion. [This last 
parameter should not be confused with the metric 7 a &, the Eddington parameter 7, nor its 
body-dependent generalization 7^4^.] Within a given theory of gravity, these three timing 
observables^ can be predicted in terms of the masses and of the pulsar and its 
companion, which are a priori unknown. The theory is consistent with experiment if there 
exists a pair of masses (mi,mg) giving the correct observed values for the three quantities 
P, Cj and 7timin g - This is the so called "P-cj-7" test, that general relativity passes with 
flying colors, and which establishes the reality of gravitational waves. We give below the 
experimental values quoted in |43|| : 



P = 27906.9807804(6) s , (5.12a) 

e = 0.6171308(4) , (5.12b) 

P = -2.422(6) x 10~ 12 , (5.12c) 

u = 4.226621(11) yr" 1 , (5.12d) 

7timin g = 4.295(2) x 10- 3 s , (5.12e) 

where figures in parentheses represent la uncertainties in the last quoted digits. In fact 
the determination of P is so precise that it is necessary to take into account the small 
variable Doppler effect due to the acceleration of the binary pulsar towards the center of 
the Galaxy ||44j| . This induces an extra contribution to P, which takes the value P ga i = 
—0.0124(64) x 10~ 12 in general relativity. The intrinsic variation of the orbital period (due 
to gravitational radiation damping) is thus given by 



18 In the language of J32LE2|, we perform a combined theory-dependent analysis of several inde- 
pendent pulsar data. 

19 In order not to create any confusion with our use of the word post-Keplerian in the present 
paper, we refer to these quantities as "timing observables" instead of "post-Keplerian" parameters 



as used in [31,32|. 
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Observed - Pgal = -2.4101(85) X 10~ 12 . (5.13) 



In tensor-scalar theories of gravity, P ga \ is modified by a small contribution, SP ga \, that we 
will take into account in our calculations below. 

The binary pulsar PSR 1534+12 has been observed only since 1991, and the present 
experimental accuracy on its parameter P is not comparable with the one of Eq. ([5.12c ) 



above. However, this pulsar is much closer to the Earth than PSR 1913+16, and it has 
been possible to measure Co and 7timin g (with very good precision), as well as two new timing 
observables, r and s, measuring the amplitude and the shape of the Shapiro time-delay 
caused by the companion. Out of these two parameters, only s is measured with precision. 
[Note that, geometrically, s = sini is the sine of the angle between the orbit and the plane of 
the sky.] As above, the three quantities Co, 7timin g an d s can be predicted as functions of the 
masses rriA, tub within a given theory of gravity. One can therefore test if this theory agrees 
with experiment by looking for a pair of masses (mi^) consistent with the three observed 
values of Co, 7timing and s ("u-j-s" test). We shall make use of the latest experimental data 
discussed in [H51: 



P = 36351.70267(2) s , (5.14a) 

e = 0.2736771(4) , (5.14b) 

x = 3.729458(2) s , (5.14c) 

Co = 1.75573(4) ° yr" 1 (5.14d) 

7timin g = 2.081(16) x 10" 3 , (5.14e) 

s = 0.981(8) . (5.14f) 

Here x = a\s/c is the projection of the semi-major axis (ai) of the pulsar orbit on the line 
of sight (in light-seconds). A precision should be given concerning the quoted experimental 
uncertainties. Ref. |45| gives (for the more reliable 1.4 GHz data) the Ax 2 = 2.30 and 
Ax 2 = 6.17 contours in the (r, c) plane with c = y/1 — s 2 . We deduced from these the 
Ax 2 = 4 contour which defines, when projected onto the c 2<7 sta t interval for c 

considered alone. We use the corresponding 2cr sta t interval for s = yl realistic 
lcr interval (in other words, we double the statistical la interval for s to take into account 
possible systematic effects). We similarly doubled the statistical la uncertainties obtained 
in ||5J for Co and 7timin g - 

The binary pulsar PSR 0655+64 is composed of a neutron star of mass ~ 1.4m , and a 
light companion of mass ~ O.8m , which is probably a white dwarf. The gravitational waves 
emitted by such a dissymmetrical system involve a large dipolar contribution of order 1/c 3 in 
tensor-scalar theories, whereas the dominant radiation in general relativity is quadrupolar 
and of order 1/c 5 . The fact that the observed value of P is very small (and consistent with 
zero) constrains therefore the existence of scalar fields, or more precisely the magnitude of 
their interaction with matter. We will see below that this system imposes a tight bound on 
the 2PN parameter The experimental data that we will need for our analysis are taken 



from 45 



P = 88877.06194(4) s , (5.15a) 
e < 3 x 10~ 5 , (5.15b) 
P = (1 ± 4) x 10" 13 . (5.15c) 
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The masses of the pulsar and its companion are not known independently; several 
pairs (mA,m B ) are thus a priori possible, such as (1.3Om , O.7m ), (1.35m , 0.8m ) or 
(1.4Om , O.9m ). In our calculations below, we will choose the mass pair which gives the 
most conservative bounds on the 2PN parameters, namely itia = 1.3Om , m B = O.7m . 
[Smaller masses could be consistent with experimental data, but current lore favors neutron 
star masses close to 1.4m .] 

The fourth and last binary pulsar that we will consider, PSR 1800—27, is also a dis- 
symmetrical system, involving a neutron star of mass rriA ~ 1.4m Q and a light companion 
of negligible self-energy. The acceleration of the pulsar towards the center of the Galaxy is 
therefore proportional to the self-gravity- modified effective gravitational constant Gao [cf. 
Eqs. (|4.20a|) and (|4.23a|) above], whereas the companion is accelerated by a force propor- 
tional to the weak- field gravitational constant Goo — G, Eq. Q2.14 ). As shown in [46], this 



violation of the strong equivalence principle causes a "gravitational Stark effect" on the orbit 
of the system, polarizing its eccentricity in a particular direction. A highly circular orbit, 
like the one of PSR 1800—27, is therefore very improbable. A statistical study can thus be 
performed to constrain the magnitude of the matter-scalar field interaction. Following the 
method of J46|, Ref. J45| has obtained the bound 



\5 A \ < 1.4 x 10~ 3 (5.16) 

on the parameter 5a characterizing this gravitational Stark effect. This bound corresponds 
to the 90% confidence level, which plays the role of an "effective la level" for the non- 
Gaussian statistics of this test. [Twice this value gives the 95% C.L., i.e., the standard 2a 
level.] Note that it is more secure than the standard la (68%) level. 

The analytic expressions of all the observable quantities discussed above have been de- 
rived in |46lj32]j24| . The theoretical prediction for the observed time derivative of the orbital 
period has the form 



pmonopolc . pdipole pquadrupolc . pquadrupolc , p , rn _|_ /O / X 1 fKT7\ 
"spinO + spinO + ^ spinO + * spin2 ' + * gal + O-Tgal + U I — | , ^0.± /J 



where the different contributions are given in Eqs. (6.52) and (9.22) of f24" |. The same refer- 
ence gives also the expressions of the periastron shift Cj and of the time-dilation parameter 
7timin g in Eqs. (9.20), as well as the "Stark" parameter 5a in Eq. (9.16). Finally, the theo- 
retical prediction for the timing observable s is given in Eq. (3.15) of [0^| . Introducing the 
notation 

M = m A + m B , X A = m A /M , X B = m B /M , n = 2n/P, (5.18) 
it can be written as 

- = — (G AB Mn/c 3 ) 1/3 ■ (5.19) 
s n 

These theoretical predictions give (when written out in detail) the various timing observables 
as theory-dependent functions of the masses of the bodies. In general, because of possible 
nonperturbative strong-field effects |4(J, the latter functions should be considered as junc- 



tionals of the 1 + n(n — l)/2 arbitrary functions entering the definition of a tensor-scalar 
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theory (notably A(<p)). However, if we assume the absence of genuine nonperturbative ef- 
fects, we can expand the functions giving the observables in powers of the compactnesses c A , 
c B of the bodies. This has the effect of reducing the functional dependence of the observables 
to a dependence upon a finite number of theory parameters. At the lowest orders, there 
appear the 1PN Eddington parameters (3 and 7. As these are already tightly constrained 
by solar-system data, Eqs. (|5.1| ), we will neglect them in the following, and investigate the 
limits that can be set on further theory parameters. The 2PN parameters e and ( appear 
precisely at the next significant order of the expansion in compactnesses pi |. The deeper 
layers of theory-parameters introduced in G3J always appear with higher powers of the com- 



pactnesses. [As shown in Table 3 of |24j], this is true even for the different contributions of 
the gravitational radiation, in spite of their rather complicated structure.] 

Here we shall estimate what constraints are imposed by binary-pulsar data on e and £, 
when restricting our attention to the lowest-order terms in compactnesses involving them. 
From a numerical point of view, our approximation is rather rough, as the next orders that 
we neglect are only 0(ca) ~ 0.3 smaller than the terms we will consider. A way of justifying 
our approach is to say that we assume that there is no fine-tuned compensation between 
the e and ^-dependent terms we retain and the higher-order terms which involve new (3PN, 
4PN, . . . ) parameters. We believe that our simplified analysis will give at least the right 
order of magnitude for the constraints on e and £ that would follow from a more complete 
theory-dependent analysis. 

To write in closed form the truncated expressions of the different observables, let us 
introduce the notation 20 



V = (G*Mn) 1/3 



(5.20) 



where M and n have been defined in Eq. (|5.18|) above, and where G* is the bare gravitational 
constant. Note that since we neglect here a 2 , oc 7, the bare constant can be identified 
with the gravitational constant G = G*(l + a 2 ,) measured in weak-field conditions. The 
timing observables can now be written as 



pmonopole 
spinO 



e 2y/2 



(5.21a) 
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1 + e 2 /2 
1 - e 2 ) 5 / 2 



+ 0(s E 



{X A -X B )(c A -c B )( 



1 + 3e 2 + 3e 4 /8 
(1 - e 2 ) 7 / 2 



32tt V £ 



X A X B [(c B X A + c A X B )\ + (b B X A + b A X B )(e + 2() + 0(s' 



1 + 73e 2 /24 + 37e 4 /96 



[1 - e 2 ) 7 / 2 



(5.21b) 



(5.21c) 



192vr V £ 
5 ~& 



X A X 



B 



1 + \c A c B Q + ^(b A + 6 B ) (e + 2C) + 0(s 3 ) 



20 The notation V is a reminder of the fact that this quantity measures some mean orbital velocity. 
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(5.21d) 

(5.21e) 
(5.21f) 

(5.21g) 
(5.21h) 

(5.21i) 



As in Eqs. ( 3.27|) above, s in the error terms on the right-hand side is a global notation for the 
compactnesses of the bodies, which should not be confused with the timing observable s on 
the left-hand side of (|5.21i| ). The 0(V 5 /c 5 ) contribution in P^^ e can °f course be neglected 
with respect to its 0(V 3 /c 3 ) term. On the other hand, the monopolar and quadrupolar 
contributions should not be neglected in spite of their being also of order 0(V 5 /c 5 ). Indeed, 
the dipolar term ( |5.21b| ) can become very small if the pulsar A and its companion B are 
almost identical, and the monopolar and quadrupolar contribution then dominate. The 
galactic contribution to P given in Eq. ( |5.21e| ) corresponds to the case of PSR 1913+16, and 
should not be used for other pulsars. The parameter k a entering (|5.21g|) is the logarithmic 
derivative of the inertia moment I A of the pulsar with respect to the gravitational constant: 
k a = d\nI A / d\nG. It has been estimated in Ref. ||7| to range between 0.5 and 1.7, 
depending on the nuclear equation of state used to describe the neutron star matter. The 
"P-cj-7" test in PSR 1913+16 is almost insensitive^] to the value chosen for k a . On the 
contrary, the "c+7-s" test in PSR 1534+12 gives slightly tighter bounds on e and ( for 
k a = 1.7 than for k a = 0.5. In order to derive conservative bounds for these parameters, 
we will use the value k a = 0.5 in the following. 

the compactnesses of the bodies can be estimated by using 
For a realistic equation of state of matter inside a neutron 



In all the equations ( |5.21 



the results of Appendix B of p4 
star, we found in this reference 



c A ps 0.21 m A /m Q , 
a A ~ 2.16 c 2 A , 
b A ~ 1.03 c\ . 



(5.22a) 
(5.22b) 
(5.22c) 



The coefficient 0.21 of Eq. (|5.22a|) can be lowered to ~ 0.15 for a stiff equation of state, and 
increased to ~ 0.30 for a soft equation of state. We will choose the central value 0.21 in our 
calculations. When the companion of the pulsar is not itself a neutron star, but a weakly 
self-gravitating body like a white dwarf, its compactness c B c A can be neglected. This is 
the case for PSR 0655+64 and PSR 1800-27. 



21 



This follows from the near parallelism of the P and ui curves in the mass plane. 
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The predictions (|5.21|) of tensor-scalar theories can now be confronted to the experimen- 
tal data Q5.121 ) — (|5.16|) . Since we are working at the first order in e and (, we can replace 



the masses tua, m B of the bodies (as well as the corresponding compactnesses ca, cb, a a, 
ag, bA, bp) by their general relativistic predictions in all the terms involving one of these 
parameters. By contrast, the dominant contributions in ^,^2 rup ° e , 7timin g and x/s should 
be considered as functions of two a priori unknown masses tua, n^B- For each value of e and 
C, we can thus determine if the above four tests can be passed. More precisely, in the cases of 
PSR 0655+64 or PSR 1800—27, the measurement of the observables P or 5a directly defines 
(when assuming the above indicated values of the masses) a one-sigma constraint on some 
linear combination of ( and e. In the cases of PSR 1913+16 or PSR 1534+12, the simulta- 
neous measurement of three observables (P-c+7 or w-j-s), which are predicted to be some 
functions of the four quantities rriA, tub, C an d £, defines (by eliminating the two unknown 
masses between the three equations and by adding in quadrature the la errors on the three 
observables) a one-sigma constraint^ on some linear combination of ( and e. Summarizing: 
each set of pulsar data leads to a reduced x 2 °f the form Xp((, £ ) — (C + ^p £ ~ f i p) 2 / cr p: 
equivalent to the one-sigma constraint — op < ( + Xpe — jip < op. We find the following 
bounds (at the la level for the first three tests, and at the 90% CL. for PSR 1800—27) 



PSR 1913+16 
PSR 1534+12 
PSR 0655+64 
PSR 1800-27 



-4 x 10" 4 < C - 5 x 10- 2 e < 7 x 10~ 3 , (5.23a) 

-8 x 10~ 2 < C + 0.15e < -10~ 4 , (5.23b) 

-7 x 10~ 3 < C < 4 x 10" 3 , (5.23c) 

|C + e/2| < 1.5 x 10~ 2 . (5.23d) 



These four allowed regions of the e-£ plane are displayed in Fig. || Clearly pulsar data favor 
only a small neighborhood of the origin e — £ = 0, i.e., of general relativity. To combine 
the constraints on e and ( coming from different pulsar experiments, we have added their 
individual \ 2 ( as defined above) as if they were part of a total experiment with uncorrelated 

Gaussian errors: xL^C) = X?9i3+i60>0 + X? 5 34+i2^ 5 C) + Xo655+64^>0 + X?soo-270> 0- 
[In spite of the non-Gaussian statistics of the gravitational Stark effect in PSR 1800—27, 
we used the bound ( |5.16| ) as an "effective la level".] In the approximation explained above, 
each x 2 is quadratic in e and (. Therefore, the sum Xtot( £ > C) is a quadratic form in e and £. 
The contour level Ax t 2 ot (^> C) = 2.3, whereQ Ax t 2 ot 

— Xtot (Xtot)mm; defines for two degrees 
of freedom the 68% C.L. (la level) ellipse represented in Fig. |9|. 



22 An alternative method is to start from the full x 2 ( m A, rnp, C, e) = ^aiQ^ ~ Qa( m A,'mB, 
C, e)] 2 /a" 2 associated with the three measurements q a = g° bs ± a a , and to reduce it to a function of 
C and e by minimizing over wia and mg. Note that we neglect the correlations between the three 
observables q a . 

23 The nice global consistency of the independent pulsar tests, proven by the overlapping of the 
strips in Fig. || means that the overall goodness-of-fit criterion associated with (xtot)min is satisfied. 
This entitles us to use the variation of xiot above (x 2 t)min to define meaningful error levels on e 
and C- 
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In conclusion, our analysis of these four binary-pulsar tests yields the bounds (at the 
combined 68% C.L.) 

\e\<7xl0~ 2 , |C| < 6 x 10 3 . (5.24) 

Because of our truncation of the observables ( |5.21| ) to their lowest order term in e and (, these 
values should be considered only as estimates of the constraints that binary-pulsar data can 
provide. They show nevertheless that possible 2PN deviations from general relativity can 
be tested with great accuracy in binary-pulsar experiments, whereas we saw in subsection 
V-B that they are almost impossible to detect in the solar system. 



VI. CONCLUSIONS 



We proposed in this paper a new, theory-based framework for conceiving and interpret- 
ing experimental tests of relativistic gravity. Previous frameworks were characterized by 
a phenomenological attitude. Eddington [JXJ] initiated such an approach by assuming that 
the (static) spherically symmetric one-body solution in a generalized relativistic theory of 
gravity could differ from the Schwarzschild solution in having arbitrary coefficients in front 
of the different powers of the small parameter Gm/pc 2 . Namely, he wrote 



#00 

9ij 



1 



Gm „ ( Gm 
2a — + 2(3 



pc* 



pc z 



+ 



Gm 



l + 2 7 — + 

pc z 



(6.1a) 
(6.1b) 



thereby introducing at the first post-Newtonian (1PN) level two independent phenomeno- 
logical parameters f3 and 7 with values one in general relativity [after having remarked that 
the Newtonian level parameter a can be conventionally fixed to unity]. The Eddington ap- 
proach has been extended in several directions. Following some pioneering work of Schiff |fj| 
and Baierlein [fl, Nordtvedt || and Will || introduced ten independent phenomenological 
parameters, j3, 7, £, ax, a 2 , a 3 , Ci> C2, C3; C4; to describe the most general N-body met- 
ric at the 1PN level. [However, their subsequent work made it clear that (3 and 7 played 
a key role among the ten parameters of their extended "PPN" formalism.] Epstein and 
Shapiro [|12|], and Fischbach and Freeman ||13| , extended the original Eddington expansion 
( |6.1|) of the one-body metric by introducing a parameter €es = &ff describing the second 
post- Newtonian (2PN) order contribution to the (isotropic) spatial metric: 



9ij — o%j 



l + 2 7 



Gm 



pc 



2 +2 e ES 



Gm 
pc 2 



+ 



(6.2) 



Benacquista and Nordtvedt [f20|-p2]1 tried to extend directly the N-body PPN formalism to 
the second post-Newtonian order by introducing a large number of a priori independent 
parameters. Finally, in a somewhat different vein, Damour and Taylor |31~l,|32f introduced a 



phenomenological approach specifically adapted to extracting the maximum possible number 
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of relativistic gravity tests from binary pulsar data ( "parametrized post-Keplerian" formal- 
ism). 

By contrast with such phenomenological approaches, we have systematically adopted in 
the present paper (which is an extension of our previous work |[24|| ) a theory-based approach. 



Instead of trying to parametrize any conceivable phenomenological deviation from general 
relativity, we work within the simplest and best motivated class of non-Einsteinian theories: 
the tensor-multi-scalar theories in which gravity is mediated by a tensor field {g^ u ) together 
with one or several scalar fields (ip a ; a = 1,2, ... ,n). These theories are, in our opinion, 
preferred for three types of reasons: (i) massless scalars naturally appear as partners of 
the graviton in most unified theories (from Kaluza-Klein to string theory); (ii) this is the 
only known class of theories respecting the basic tenets of field theory (notably the absence 
of negative-energy excitations) in which very high precision tests of the equivalence prin- 
ciple can naturally be compatible with post-Newtonian deviations at a measurable level^; 
(iii) they naturally "explain" the key role played by the original Eddington parameters f3 and 
7 at the 1PN level, and have a simple enough structure (in spite of their great generality^) 
to allow one to work out in detail their observational consequences at the 2PN level. 

Our main results are the following. Two, and only two, new parameters (beyond the 
non-Einsteinian 1PN parameters /? = /? — 1,7 = 7 — 1) quantify possible non-Einsteinian 
effects at the 2PN level: e and (. The role of these 2PN parameters is threefold: 

(i) They parametrize 2PN deviations from the general relativistic (iV-body) physical 
metric tensor: 

+ | !W /«^ + o(|I) + o(I), 

5 goi (x) = o(^,l) +o(Jf) , (6.3b) 

^) = °(?'?) + °Qf) ■ ( 6 ' 3c ) 

(ii) They determine the renormalizations, due to self-gravity effects, of the various effec- 
tive gravitational coupling parameters between massive bodies: 



24 For instance, in all (positive-energy) vector theories coupled (at the linear level) to some current, 
the existing equivalence principle tests at the 10~ 12 level necessarily constrain all post-Newtonian 
deviations at the ~ 10 -9 level. Though current unified models suggest the existence of massless 
scalar fields with composition-dependent couplings, they leave open the possibility of a very small 
parameter (e.g. a factor ~ 1CP 5 was found in string-inspired models p6[| ) relating post-Newtonian 
deviations to equivalence-principle tests. 

25 Indeed, we consider the most general n-scalar models described by 1 + n(n — l)/2 arbitrary 
functions of n variables. 
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\rriAC z itlbC a J \rriAC z J \m,BC J 

+ (| + C) {U2)A + {U2)b + 0{M) x 0(s 2 ) + 0(, 3 ) , (6.4a) 
Jab = 7 + O^, 7) x 0(5) + 0(s 2 ) , (6.4b) 

/? fl c = + (e + + C H4 + -S + °(P> 7) x + °(* 2 ) • ( 6 - 4c ) 

m^cr \rriBC z mcc z J 

Here, Gab denotes the effective Newtonian constant measuring the strength of the 

A 

0{vn avtib I tab) leading gravitational coupling between A and B, while ^f AB and j3 BC denote 
effective Eddington parameters measuring the strength of the 0{mATriB{^A — v b) 2 '/ h 'abc 2 ) 
and 0(mATnBmc / '(r ab^ acc 2 )) post-Newtonian couplings. Moreover, 77 stands for the com- 
bination 4/3 — 7 = 4/3 — 7 — 3, while s ~ E era,v /mc 2 denotes the strength of self-gravity 
effects. The more complete version of our results on self-gravity renormalizations is given in 
Eqs. flpg ). 



(iii) They determine the renormalizations of the locally measured coupling parameters 
(e.g. the local gravitational constant in physical units G\ oc ) due to the presence of distant 
"spectator" matter (say, a mass at a distance D) around the considered gravitating 
system: 



G\ oc — G c 



Gm s 

1 — T] 

Dc 2 1 



+ 0(-U , (6.5a) 



D 



7ioc = 7oo + 4^/3(2 + 7) + 0^) , (6.5b) 



Aoc = /3oc - ^ (I + C - 8/3 2 ) + O (i) . (6.5c) 

On the other hand, we find that £ and £ tio not enter light-deflection nor time-delay 
experiments at the 2PN level (see Appendix C). In particular, we find that the "post- 
Eddington" formal 2PN parameter introduced in Refs. |T2|Jl~3[| , as recalled in Eq. ( |6.2|) above, 
is the following function of the 1PN parameters (3 = (3 — 1, 7 = 7 — 1: 

e ES = 1 + ^+^7- (6.6) 

6 6 

This shows that second-order light-deflection (or time-delay) experiments do not probe any 
theoretically-motivated 2PN deviations from general relativity. More generally, after dis- 
cussing observable effects linked to e 7^ or ( 7^ (in planetary perihelia, Lunar Laser 
Ranging or other strong-equivalence-principle tests, and clock experiments), we conclude 
that the solar system is not an appropriate testing ground for probing possible 2PN devi- 
ations from general relativity. More precisely, we find that in the best cases, foreseeable 
technology might barely be able to constrain e and ( at the order unity level, while, in the 
worst cases, these parameters might complicate the interpretation of 1PN experiments by 
contaminating some observables. 

This seemingly negative conclusion is, however, to be tempered by the following two 
positive conclusions of ours: 
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(a) We identified binary pulsar experiments as an excellent testing ground for the 2PN 
structure of relativistic gravity. By a simplified (linearized) analysis of existing data on the 
binary systems PSR 1913+16, PSR 1534+12, PSR 0655+64 and PSR 1800-27, we were in 
fact able to constrain already e and ( at the level 

|e| < 7 x 10^ 2 , |C| < 6 x 10 3 . (6.7) 

(b) We stressed that solar-system experiments are well suited to measuring with high 
precision the 1PN parameter 7 which is of greatest significance among all post-Einstein pa- 
rameters. Indeed, our theory-based analysis shows very clearly that 7 is a direct measure of 
the coupling strength of matter to the scalar fields. All the other post-Einstein parameters 
(P, £,£,...) necessarily tend to zero with 7 in all theories having only positive-energy exci- 
tations (more about this below). From this (theory-based) point of view, the most important 
solar-system experiments would be high-precision light-deflection or time-delay experiments 
reaching the level, say, 7 < 3 x 10 -7 , which comes out naturally from mechanisms of cos- 
mo logical attraction of tensor-scalar models toward general relativity |p7|j26l 26 . 



Let us end by stressing again some of the conceptual and technical differences between 
the theory-based framework introduced here, and the usual Eddington-Nordtvedt-Will PPN 
framework (and its extensions). The natural range of values of the PPN parameters is 
uniformly supposed to be of order unity, independently of the post-Newtonian order at 
which they appear. In our approach, there is one basic set of (experimentally small) coupling 
parameters, 

where m(ip) is the mass of a particle in the Einstein conformal frame. [This way of writing 
the definition of a a is of great generality as it encompasses both self-gravity effects |24 and 
possible composition-dependent effects [p^] .] 

All the phenomenological parameters 7, (3, (, e, . . . measuring the 1PN, 2PN, etc. de- 



viations from general relativity (see |24| for a list of the 3PN and higher parameters) are 
explicitly constructed by contracting at least two of the basic a a 's with objects built from 
successive covariant field derivatives of the a a 's: (3 a b = D a ab, /3' abc = D a Dba c , etc. (see 
section II B for the definition of D a ). In particular, internal indices in the scalar-field space 
being raised and lowered by means of the positive-definite metric 7 a b(y? c ) defining the kinetic 
energy of the scalar fields, 



a n a a 



7 = -2 - a , (6.9a) 

' l + a a a a v ; 

- 1 a a p ab a b 

P = 2 (l + a a «»)* ' (6 ' 9b) 



26 We note that the space missions GAIA (Global Astrometric Interferometer for Astrophysics) 
and SORT (Solar Orbit Relativity Test), proposed by the European Space Agency, aim at such a 
level for 7. 
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(1 + a a a a f ' 
3\ a a a b a c 

e = ^ abc " . (6.9d) 
(l + a a a a ) 3 v ' 

Therefore, contrary to the PPN philosophy where, say, the 1PN parameters 7 and (3 could 
be accidentally small and the 2PN ones ( and e of order unity, our approach suggests that 
all of them tend to zero with 7, or more precisely with a 2 = —7/ (2 + 7) which is a positive 
measure of the total coupling strength of matter to the scalar field. [The ones, like e, which 
are cubic in the a's are even expected to be much smaller than the quadratically small ones^] 
7, P, C-] 

Finally, let us note that, as far as we are aware, the present work is the first one to 
use in an effective way a direct classical diagrammatic approach to the relativistic iV-body 
problem^ We hope to have convinced the reader of the technical power of this method. 
Indeed, once one is used to the notation, our results ( |6.9|) on the complete set of parameters 
entering the first two post-Newtonian levels can be obtained in drawing just half a page of 
simple diagrams. 
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APPENDIX A: EXPLICIT DIAGRAMMATIC CALCULATIONS 

In order to compute explicitly the diagrams of Fig. [7], we must first derive the expressions 
of the propagators and of the different vertices of Fig. f|. As discussed in subsection III-A, we 
need to expand the gauged- fixed action of the theory up to the fourth order in = g^ — f^v 
and ip a — ifQ. We need only the terms quadratic in in the Einstein-Hilbert action fl2.4p 
to define the graviton propagator (see e.g. for higher-order terms). We get easily 



c 3 



*Sspin2 — 1 „ ^ I d X 
l07TG* 



(Al) 



where the indices are raised with the flat metric = diag(— 1, 1, 1, 1), and where Q al3 "' 6 = 
\{f ai f 136 + f a6 f^ 1 ~ / a/3 / 7<5 )- To invert this kinetic term, it is necessary to fix the gauge. 



27 Note that the expected extra smallness of e gives a mnemonic rule for remembering its definition. 
28 Bertotti and Plebanski [37| mainly discussed general features of a classical diagrammatic ex- 



pansion, while previous work by the Japanese school |ljj,[48| had used quantum diagrams and then 



converted S-matrix elements into an effective iV-body potential. 



46 



The most convenient choice for our calculations will be the harmonic gauge, corresponding 
to a gauge-fixing term — (c 3 /327rG*) f d A x(d v h u ^ — |<9 M /i^) 2 + 0{h 3 ) (see e.g. Jl7| for the exact 
harmonic-gauge-fixing term). This cancels exactly the second term of (|A1), and we get after 
an integration by parts 



S&a = J d A x l -h aP Q^< s U f h l8 + tot.div. + 0(h 3 ) , (A2) 

where □/ = j^ v Qpd v is the flat d'Alembertian. This kinetic term has therefore the form 
f d A x(—^hV^h), where Vh defines the graviton propagator, represented as a curly line in 
Fig. [7|. In terms of the Green function fl3.17| ), satisfying DfQ(x) = — An5^(x), we thus get 

V^ s (x, y) = -^-PaMS Qi? ~ V) , (A3) 

where P aPl& = f aJ fp S + fasffr ~ fapf~ ( 5 is the inverse of Q a ^ s . 

Let us now expand the action ( |2.5| ) of the scalar fields up to the fourth order in (p a — y>g 
or hfj, v . To simplify the expressions, it is convenient to choose Riemann normal coordinates 
at ip in the internal scalar space, so that the metric 7 a b can be expanded as in Eq. ( |3.15 ). 
It is also useful to define the origin of this scalar space at (p , i.e., to choose the coordinates 
ip a of this space so that y?o = 0. Then, the scalar-field action ( |2.5|) can be expanded as 



Sspino = J d A x \d^ a d v y h (r - hT + W p h p » + 0(h 3 )) 

X (jabfao) - ~#ac6d(^o)^V + 0(ip 3 )^j , (A4) 

where the terms inside parentheses are respectively the expansions of g$ v , Jg~l, and r y a b( l p)- 
The kinetic term of the scalar fields reads therefore (c 3 /47rG*) J d A x ^tp a 'y a b(tp )\-\f(p b + 
tot.div. = J d A x{—\ipV~ l ip) ) where 



K h &y) = ^ ab MG(x-y) (as) 



c 



is the scalar propagator, represented as a straight line in Fig. [7[ [As before, j ab denotes the 
inverse of 7^.] 

We can also derive from (A"4) the expressions of the vertices connecting scalar fields and 
gravitons. Our conventions for defining vertices are the following: (i) We first define some 
formal "global" vertices Vi = iSi when considering (as in Figs. and |j) the gravitational 
sector as a whole, $ = (</?, h). (ii) We then define the individual vertices Vi(ip, h) by formally 
expanding the global multilinear forms V^($) as if <£> were equal to the sum <p + h, e.g. 

V 3 (<p + h,(p + h,(p + h) = V 3 ((p, <p, <p) + W 3 ((p, <p, h) + 3V 3 (<p, h, h) + V 3 (h, h, h) . (A6) 

This convention allows us to use directly the simple global diagrams of Fig. [] with the 
intuitive replacements of Fig. [7]. For instance, the vertex of order 0(ipiph), entering the first 
T diagram of Fig. [7| as well as the first four F diagrams and the first three H diagrams, reads 
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T wh = — J d 4 x jaMQ^d^d^Kp , (A7) 

where Q^ va ^ is the same tensor as in ( |A1| ) above. In other words, the kernel defining the 
T w /, vertex is 

3 

(T w ,)f (x 1; x 2 ,x 3 ) = -^G;^o)Q^g^6(x 3 - xi)-^S(x 3 - x 2 ) . (A8) 

Note that it is symmetric in x\, x 2 (scalar lines), but not in all variables since they represent 
physically inequivalent lines. The vertex connecting two scalar lines and two graviton lines, 
entering the second X diagram of Fig. 0, reads 

Xwhh = ~t^g~ I dA%lab ^ {^r a Q vPlS - l -rQ a ^ d^ a d vV b K P h l5 . (A9) 

Finally, the vertex connecting four scalar lines, entering the first X diagram of Fig. [I], reads 
X„ 4 = J d 4 xR acbd (v )^%<p a d^ b , (AlO) 



where R a bcd is the Riemann curvature tensor of jab- It should be noted that the four scalar 
lines of this vertex are not equivalent, since two of them involve the derivatives d li (p a of the 
fields. The same remark holds also for the multi-graviton vertices of order 0(h 3 ) and 0(h A ) 
in ([A~l]). One should therefore symmetrize these vertices [i.e., write the distributional kernel 
X,p4,(xx, X2, X3, X4) read off (|A10|) as a symmetric function of X\, x 2 , £3, x&] before computing 
the Fokker action ( |3.14j ). Alternatively, one can also use their non-symmetric form, but take 
into account the different ways to choose the lines involving derivatives in the diagrams of 
Fig. |7|. This does not change anything for the X diagrams, since the Lagrangian is anyway 
symmetrized over the 4 bodies A, B, C, D, but this leads to the numerical weights displayed 
in Fig. |H] for the F and H diagrams. 

Let us now consider the matter action (|2.16|) , describing N self-gravitating bodies. As 
shown in subsection III-A, we need to expand it only up to the third order in h^ u and (p a . 
Still using Riemann normal coordinates at (p = in the </>space, we get easily 



AT , 



S m = I m° A c 2 dT A 



A=l 
1 



+ g(fl& + P>c + PL< + #X + a*a£a? ) <A>V + <% 4 ) 



(All) 



where the first bracket comes from the expansion of m^tp) around <y9 , and the second one 
from the expansion of ds* A around the flat metric. All the fields appearing in ( |A11| ) must 
be evaluated on the worldline x M = x^(ta). As above, dr A = (1 — v^/c 2 ) 1 / 2 ^ denotes the 
Minkowski proper time of body A, u A = dx A /cdrA denotes its (Minkowski) unit 4- velocity, 
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and the index o means that the corresponding quantities are evaluated at (p a = <^g (= 0). The 
different vertices connecting material bodies to gravitons or scalar fields can thus be read 
directly from ([All ), taking into account both the factors 1, 2 or 3 entering their definition 



in Fig. f| and the binomial factors coming from our convention (ii) above, cf. Eq. (|A6|). The 
linear interaction terms read 

l v = -J2 I dT A m° A c 2 (a^) ip a (x A ) , (A12a) 

A J 

h = / dr A m A c 2 ]-u A ^ l3 A h a p(x A ) ■ (Al2b) 

A 2 

The corresponding spacetime sources (white blobs) read thus explicitly 

[ dr A m A c 2 {a%8^{x - x A {r)) , (A13a) 

A J 

\Y,J dr A m A c 2 u a A u^ A \x - x A (r)) . (A13b) 



(TAX) 



a a/3 (x) } 

~ A 



As for the nonlinear interaction with the source (vertices connecting matter to several field 
lines), they read, when omitting a common J2 A I d T A m A c 2 5^(x — x a (t)) in front: 

(A14a) 
(Al4b) 

(Al4c) 

, (M4d) 
(A14e) 
(Al4f) 
(Al4g) 

The calculation of the different diagrams of Fig. [7] can now be performed straightfor- 
wardly. Let us summarize our diagrammatic rules: The Fokker action is a sum of contri- 
butions, each of which is represented by a diagram endowed with a numerical coefficient 
obtained by combining the factors indicated in Fig. ^ and in Fig. [7|. Each (bare) diagram is 
computed by the following rules: (i) replace the white blobs by Eqs. ( |A13| ) if they have only 
one incident line, or by Eqs. ( |A14| ) [with an extra ^2 A J dT A m A c 2 8^\x — x A {r))} if they have 
several incident lines; (ii) replace each internal line by the appropriate propagator V(x,y); 
(iii) replace each field vertex by a suitably symmetrized distributional kernel; (iv) integrate 
over all the spacetime points. An additional rule is that infinite self-interactions are dis- 
carded. (As they have the same structure as in general relativity, this rule can probably be 
justified by methods similar to the ones used in deriving the 2PN Lagrangian in Einstein's 



V w : 






2 a a U >A , 




1 a P 7 S 

-u A u A u A u A , 








\{^ b + atat)u a A u A 




\^u A u A u\u A , 


£hhh '■ 


3 a /3 7 8 e f 

—u A u A u A u A u A u A . 
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theory |14j| .) In order to save space, we will not compute each of the diagrams in the present 
appendix, but rather show the technique on the most important of them. 
Let us start with the 2-body interaction term 1 1 of Eq. ( |3.14| ). We get 

\^=\ff dxdy [a a (x)V;\x,y)a b (y) + a^(x)r^ s (x,y)a^(y)] 

= \ ^2 [ dTA [ dl ~ B ( m W)( m B c2 )-rG( x A - x B ) [{a A a B )o + 2(u A u B ) 2 - l] , (A15) 

where (a A as) = ^al ah<x b comes from the scalar propagator (first I diagram of Fig. 0), and 
2{uaUb) 2 — 1 = u^u 1 ^ Pa^s u B u 5 B comes from the graviton propagator (second I diagram). 
Introducing the notation Gab = G*[l + (0:^05)0] and 2 + j AB = 2/[l + (0^0:5)0] as in 
Eqs. ( |2.21| ) above, we thus get 

-1 = ~y^[ dr A I dr B G AB m\m B [l + (2 + 7ab) {{u A u B ) 2 - l)] cQ{x A - x B ) . (A16) 
1 1 A+B J J 

This action describes the 2-body interaction of self-gravitating bodies having arbitrary ve- 
locities. The second post-Keplerian approximation can now be obtained by expanding this 
result in powers of v/c, using Eqs. ([3.19|) - fl3.21|) above. We get an expression of the form 
( 3.231 ), where the corrections proportional to 7ab/c 4 involve the function 

fi(r AB , va, vb, a A , a B , a B ) = [(v A - w B f (v A + v|) /2 - (v A x v B ) 2 ] /r AB 

[|x A -x B (t B )| (v A -v B (t B )) 2 ] tfl=i . (A17) 

Note that only x B and v B are time-differentiated in the second term, before setting t B = 
t. The derivative of a B can be eliminated by an integration by parts, and up to a total 
derivative, we finally get 

2 f \2 

. . Viq r 2 2 / \i \i \y A ^ V B J 
/i(r AB , v A , v B ,a A ,a B ) = [v A + v B + v A • w B - (v A • n AB )(v B • n AB )J 

^ AB r AB 

+ (v A • n AB )(a B ■ v AB ) + (v B • n AB )(a A • v AB ) + (a A • a B )r AB , 

(A18) 

where v AjB = v A - v B and n AB = (x A - x B )/r AjB . 

Let us now consider the 3-body interaction, written as |V + |T in ( p.!4|) and Fig. |6|. The 
three V diagrams of Fig. [7| give the contribution 

1 1 . f f f G 2 

-V = - dr A dr B dT C (m A c 2 )(m B c 2 )(m c c 2 )^g(x A - x B )G(x A - xc) 

b+a+c-' 

x J(2(m a m b ) 2 - 1 - (a A a B ) ) (2(u A u c ) 2 - 1 - (a A a c )o) 

-(a B p A a c )o - 2(a A a B ) (a A a c )o ■ (A19) 
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A 

Introducing as in Eqs. ( |2.21| ) the notation Gab, Iab an d Pbc-i this contribution can therefore 
be written as 

,0 ™0 ^,0 



Y r r f G a G ui in 7Tb 

-V = ^2 dr A dr B / dr c — ^ ° 2 A B c cQ(x A - x B ) cQ(x A - x c ) 

B^A^C* J ^ '' 

x { [1 + Iab + (2 + Iab) {{u A u B f - l)] [l + Iac + (2 + Iac) {(u A u c ) 2 - l)] 



-2^-2^2}. (A20) 

The 2PK approximation can now be obtained easily by expanding this expression in powers 
of v/c. Note that the Newtonian approximation of the Green function fl3.21|) , cQ{xa~ xb) — 
S(tA — ^b) / r ab + 0(l/c 2 ), is sufficient for the terms proportional to (uaUb) 2 — 1 = (va — 
vb) 2 /c 2 + 0(l/c 4 ), or (uaUc) 2 — 1. The accelerations introduced by these Green functions 
in the other terms can be eliminated by suitable integrations by parts. For instance, one 
can write 

1 Id 2 (v AB • n AB )(v c • n AC ) 

\x A - xc(tc)\ = tot.div. - 



| v A vg K V A ^ACjyv_C "-AC J (A21) 



r AB 2c 2 dt 2 c ^ 2r 2 AB c 

v A • v c - (V4 ■ n AC ) (v c ■ n AC ) 

2r A BTACC 2 

The T diagrams of Fig. [7] are a little more subtle to compute, because the central vertex 
is not located on a material body, and also because it involves derivatives of the fields. 
For instance, the first T diagram, involving one graviton and two scalar lines, gives the 
contribution 

T i = E [ dr 4 d ^B [ dr c [ d'x G " m Y B< (a A a B )o 

A,B,C^ J J J 71 

«^ 9 %;; a)9S %- xb) ^-zc), <a22) 

where x M is the arbitrary spacetime location of the vertex. The lowest order term of the 
post-Keplerian approximation reads thus 

G 2 m A m B m c 



Ti = £ / dt A J dt B J dt c J dt I d 3 x G y (a A a B ) 



A,B,C 
X U r U r t^- 

c c dx^ 



6{t - t A ) 



t x a 



_d_ 

dx L 



8{t - t B ) 



r X B 



r x c V \c 2 



A,B,C 



-Y dt d^ ^ABlABm A mBmc {-ac ; y^ac ■ n xB ) Q 1 
^ I J 47TC 4 r 2 A r 2 r x c V c 6 / 



where we have used ^^[^(t — t^)/^/!] = S(t — t A ) v ac ■ ^ x aI t2 x a c + 0(l/c 3 ). Therefore, 
although this T diagram is of the same formal order G 2 m 3 /r 2 c 2 as the V diagrams in the 
non-linearity expansion (arbitrarily large velocities), it reduces to order 0(l/c 4 ) in the post- 
Keplerian approximation (\v/c\ <C 1). This is due to the particular form of the T w /, vertex, 
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Eq. (|A7|) , which involves a specific combination of derivatives of the fields with the inverse of 
the tensor P a ^s entering the graviton propagator (|A3|) . It should be noted that the presence 
of derivatives is not sufficient to conclude that the diagram is reduced by a factor {v/c) 2 . 
For instance, the second T diagram of Fig. |7|, involving three graviton lines, does contribute 
at the first post-Keplerian order, 



dt 



G 2 m° A m B m c 



tabTac^ 



+ 



(A24) 



although the 3-graviton vertex is also of the form hdhdh. Indeed, the dominant contribution 
to this diagram is proportional to the contraction d^Qix — XA)d^Q(x — xb) = S(t — tA)8{t — 
• n X B/rl A rl B c 2 + 0(l/c 4 ), which starts at order 1/c 2 . The sum of ( |A24| ) with the 



lowest order term of ( |A20|) gives the 1PK contribution to the 3-body interaction Lagrangian, 
displayed in Eqs. (EOOcl) or ( ^25| ). 

Let us turn now to the most important 4-body interaction terms, namely the G, Z, and X 
diagrams of Fig. [7], which involve the new parameters (|3.18| ). At the second post-Keplerian 
approximation, it is sufficient to use the instantaneous Green function cQ(xa — xb) = ${tA — 
^b) 1 7™ ab in these diagrams, and one easily gets 



1 _ 1 V- f r. G abG acG AD m A m B m C m D [, , 1 e A , 2 (l3 A ,13 A ,15 A \ 
3 G "2 2^ J M TABTACTABC' I + 3^°° + 3 V BC + PbD + PcD ) 



A^(B,C,D) 

2 t _ _ _ N 1 _ __ 

+ o ^AB + 7 AC + IAD) + o \lABlAC + 1 ABl AD + IacIad) + ~1 AbI Acl AD 



+0 



2 2 ^ 

A+B+C+D 



(A25) 



dt 



GABG B cGcDm A m B m c m D 



TABTBCrCDC* 



1 + (abcd + 2[(3 a B c + (3bd 



+ 2 ( 2 + 7bc) (lAB + 1CD + lABlCD) 



+ 



(A26) 



[Note that the e contribution to the G diagram has a non-trivial normalization: G~ 
i£G 3 m 4 /r 3 c 4 , while Z ~ £G 3 m 4 /r 3 c 4 .] The contribution of the first X diagram of Fig. [7] 



2 

reads 



^ Xi = ^ Jdt A Jdt B J dt c J dt D Jdtjd 3 



Glm A m B m c m D( n ah n cd\ 

{rtacbda A OiB a C a D) 



A,B,C,D ■ 

5(t - t c ) 5(t - t D ) 

r xC 

— T 



6ttc 4 







a, 



r xD 



5(t-t A ) 



r x A 



<9" 



'5(t-t 



B 



r xB 



+ -r 




1 < J ^ gW^j XACBD + o f 1 , . (A27) 

A.BCD J J r xA r xB r xC^xD \C 



Note that although Xacbd is antisymmetric in A, C (and B,D), and although the sum is 
taken over all possible choices of A, B, C, D, this contribution does not vanish identically 
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since the integrand (n^ • n x s) / r 1 A r xB r xC r xD is not symmetric in A,C. This contribution 
vanishes nevertheless in the 2PN approximation, cf. ( p. 32 ), or if the theory involves only 



one scalar field (R a bcd(^) = 0). The contributions of the two other X diagrams of Fig. [7| are 
computed in the same way. The last one gives the usual contribution obtained in general 
relativity, and we find for the second X diagram 

^ A,B,C,D J J 'xA'xB'^C'xD \C / 

Let us end this appendix by a brief discussion of the F and H diagrams of Fig. [I], which 
are less important since the deviations from general relativity they give at 2PN order are 
proportional to 7, 7 2 or (3, and therefore already tightly constrained by present experimental 
data, Eqs. ( |5.1| ). Is is easy to see that the first four F diagrams do not contribute at the 
second post-Keplerian level, because of the presence of the T w h vertex whose graviton is 
directly coupled to a material body. Indeed, the calculation is similar to that of the first 
T diagram, Eqs. ( |A22|) - (|A23|) , and the factor {u^d^Q) 2 oc v 2 reduces these F diagrams to 
the order G 3 m A v 2 /r 3 c 6 . Similarly, the second H diagram is even of order G 3 m 4 v 4 /r 3 c 8 , and 
contributes only at the fourth post-Keplerian level. On the contrary, the first and third 
H diagrams do contribute at the 2PK level, because the graviton of their T w h vertices is 
connected to a second vertex, and not directly to a material body. These diagrams involve 



therefore contractions of the form d^Qd^Q, like in (|A24|) , and are indeed of order G 3 m 4 /r"c 



3^4 



Note that the four material bodies of these H diagrams are not supposed to be necessarily 
different from each other, since they are not directly connected by a propagator (as opposed 
to the I, V, e, or Z diagrams). Together with the first T diagram ( |A23| ) and the second X 
diagram (|A28|) , they thus yield contributions proportional to 

(rAA 7 AA = -2(a A a A ) . (A29) 



1 + j AA /2 a 

This explains the presence of this factor in the static and spherically symmetric solution 
discussed in section IV, Eq. ( }4.21|) . [The isotropic form (|4.18| ) contains extra a 2 A terms due 



to the change of coordinates ( |4.22| ).1 In fact, the contribution of the first H diagram to 
this one-body metric is proportional to (a A ao){a 2 A )\ One of the white blobs, involving the 
background value a((po) = ao, corresponds to the point where the metric is computed, and 
the three other blobs correspond to body A. Note also that the first T diagram, Eq. ( A23|) , 
does not contribute to the a A terms of g 00 , since it vanishes in the static case. Finally, let 
us mention that the last F and H diagrams are the usual contributions obtained in general 
relativity, and that the fifth F diagram is obviously proportional to the second T diagram: 

F , = E l^ x 3GgW£ + ^), (A30 ) 
where ^T^' B ' C ^ is given by (|A24|) above, symmetrized over the three bodies A, B, C. 
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APPENDIX B: 2PN RENORMALIZATIONS OF COUPLING PARAMETERS, 
THE STRONG EQUIVALENCE PRINCIPLE, AND ALL THAT 



As emphasized by Nordtvedt |p2]] , the coupling parameters (G, m, /3, 7, . . . ) of Lorentz- 
invariant gravitational theories are influenced not only by the self-gravity of the interacting 
bodies, but also by the presence of distant "spectator" matter around the system. We have 
already studied the self-gravity renormalizations in the main text. In this Appendix, we 
relate our self-gravity results to the ones of Ref. |22| and show how, within the context to 
tensor-multi- scalar theories, one can derive very easily the influence of external matter. 

Let us first mention that, from the 1PK Lagrangian ( |2.19| ), one can easily give explicit 
expressions for all the mass parameters M(G), M(7), M(/3), . . . introduced in |22j| . For the 



convenience of the reader, we give in Table | a translation of Nordtvedt's notation in terms 



of our body- dependent parameters Gab, Iab an d Pbc defined in Eqs. (|2.21| ) above. As 
before an index in one of these parameters corresponds to a non- self-gravitating body, 
so that the er-model tensors a A , [3^ of Eqs. ( [2.17[ ) should be replaced by their weak-field 

counterparts Q2.10| )- p.i2| ). For instance, j3 00 = |(a/?A«)o/ [1 + («a«)o] 2 instead of ( |2.21c[ ). 

~ a 

The 2PN renormalizations of Gab, Iab an d Pbc due t° the self-energy of the bodies have 
been obtained in Eqs. (|3.29|) above. Using them in the translation of notation of Table |, 
we recover the corresponding (less complete) results of |22 



Let us now turn to the study of the influence of external matter. The effect of a distant 
spectator body on local gravitational physics can be analyzed straightforwardly in the con- 
text of tensor-multi-scalar theories. Indeed, if this spectator S is located at a distance D 
from a local system, the local background values of the scalar fields are changed from their 
values <£>q far from the spectator to 



^0 



G*a a s m s 
Dc 2 



O 



1 

D 2 



(Bl) 



All the physical quantities / which depend on these background values are therefore renor- 
malized (compared to the value they would be measured to have at infinity) by the presence 
of the spectator as 



f f _ f G*m s ^ a df 



(B2) 



In particular, if the spectator body is supposed to have a negligible self-energy (a| = «q), 
we get for the effective gravitational constant ( [2.14| ) and the Eddington parameters ( |2.15| ) 



G 

7 



G 

7 + 4 



Gm s 

1 Tj 



Dc 2 
Gm s -5 



O 



Dc 2 

<i 

Dc 2 V2 



/3(2 + 7)+0 ( rv2 



1 

LP 



Gm s (s > — 2\ / 1 
b + )+0[- 2 



(B3a) 
(B3b) 
(B3c) 
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where 77 = 4/5 — 7 as before, and where the dimensionless ratio Gms/Dc 2 may be replaced by 
its expression in physical units Grhs/Dc 2 . Equation (|B3a|) is the well-known renormalization 
of the gravitational constant derived in H. The renormalizations of 7 and (3 have also been 



studied in Eqs. (5.15) and (3.7) of [22], but they were expressed in terms of several 2PN 



parameters, instead of the simple form (|B3b|) and (|B3c|). In particular, the analysis of 



Ref. |22| did not suffice to prove that the 2PN renormalization of 7 is in fact proportional 
to (3, and therefore that it is already constrained by the 1PN experimental bounds on 
\P\ < 6 x 10" 4 . 

Of course, the renormalizations ( |B3"1) can be generalized straightforwardly to the strong- 

~ A 

field regime, by considering the body-dependent parameters Gab^o), Iab^o), Pbc(^Po): 
and a compact spectator body (ag 7^ ag). This has been done in section 7.2 of Ref. |2"4fl , 
where we analyzed the consequences of the strong equivalence principle (SEP) in tensor- 
multi-scalar theories. This principle states that local gravitational physics is totally inde- 
pendent of the presence of spectator bodies. In other words, the renormalizations ( p3|) and 
their strong-field analogues are supposed to vanish. Let us prove that the only tensor-multi- 
scalar theories containing only positive-energy excitations (7^ > 0) and satisfying the strong 
equivalence principle are perturbatively equivalent to general relativity (to all orders). In- 



deed, we showed in p4]| , using Eqs. (p3|), that any such theory must satisfy 7 = 0, among 
other relations. From Eq. (|2.15a|) and the positivity of 7 a b, this implies «q = 0. Using now 
the diagrammatically evident^ fact that any observable deviation from general relativity 
must involve at least two factors ojq (to fill the end blobs connected to scalar propagators, 
in diagrams such as Fig. [7] or any higher-order ones), we find that all non-Einsteinian terms 
necessarily vanish. It should be noted that this result cannot be extended to the pure scalar 



theories we also considered in section 7.2 of ||24[| , in which the gravitational interaction is 



mediated by one or several scalar fields without any tensorial contribution. We showed that 
the SEP implies in that case Gab = G, Jab = 7 = —2, and /3 BC = (3 = — |, i.e., that 
the theory is equivalent to Nordstr0m's theory at the 1PK level. The assumption that the 
cr-model metric is positive does not allow us to complement this result to higher orders. 
We can also prove that some assumptions made in Ref. PI1 are inconsistent within 



the tensor-scalar framework. Indeed, Ref. |H[ assumed one could work within a class of 



theories containing no dipolar radiation and still differing from general relativity at 2PN 
order. However, we can use the fact that the dominant dipolar radiation emitted by a 
binary system is proportional to 7^(0/4 — Oi B ){a b A — a b B ). Let us first assume (as we generally 
do) that 7^ is positive definite. Then the assumption that the dipolar radiation vanishes 
implies that a a A = a a B for any body A and B, and in particular that a a A = a a if we choose a 
non-self-gravitating body B (m B = const). Using now the expansion of a a A in powers of the 



compactness of body A, that we derived in Eqs. (8.3) and (8.4) of pl| , and the fact that 
a a A = a a must be verified for any body A, we can conclude that a a = 0. In other words, 
the scalar fields are totally decoupled from matter, and the theory is strictly equivalent 
to general relativity. Even if we drop the assumption of a positive definite 7 a b, i.e., if we 
phenomenologically allow the presence of scalar fields with negative energy (ghost modes), 



29 Note that this would not be true if we were considering quantum (loop) diagrams. 
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then the same Eqs. (8.3) and (8.4) of [24] can be used to show that the theory is equivalent 
to general relativity up to the 2PN order included, but not beyond (which means that it can 
differ from it at the 2PK order, as illustrated in section 9 of ||24|| ). 

The ease with which we derived the results of this appendix illustrates again the power 
of our field-theoretical approach. 



APPENDIX C: EXPLICIT EXPRESSIONS OF THE 2PN LIGHT DEFLECTION 

AND PERIHELION SHIFT 

Before computing the explicit expression of the light-deflection angle for the one-body 



metric (|4.18| ), let us show how our diagrammatic approach allows one to get, without any 



calculations, the structure of the final result. More precisely, let us show that the light- 
deflection angle for a self-gravitating body has the structure 

O(-? AQ )(G A0 m A /rc 2 ) + O(-? A0 ^ AA )(G A0 m A /rc 2 ) 2 + 0(l/c 6 ) , (CI) 
A 

in which neither (3 BC nor the new 2PK parameters we introduced above enter. Indeed, the 
action describing the electromagnetic field minimally coupled to the physical metric is 
a conformal invariant, and can thus be written as 

Sem = -\J ^-y/Sg^F^up = -\J ^Jg-g^g^F^ . (C2) 

The second expression, involving the Einstein metric g^ u , shows that photons are coupled 
to the graviton h^ v = g*^ u — f^ u , but not to the scalar fields ip a — tp^. Therefore, the only 
diagrams describing the 2PK interaction of light with the gravitational field generated by 



material bodies are those of Fig. [11. [The G, Z, F, H and X diagrams enter the metric 
at order 0(l/c 6 ), and influence the propagation of light at the third post-Keplerian level.] 
The dominant contribution does not involve any scalar field, and is thus proportional to 
G*m A /rc 2 . However, we do not have a direct experimental access to G>, and this contribu- 
tion should be rewritten in terms of the Newtonian potential G A Qm A /r felt by a test mass 
in the vicinity of body A. This can be done thanks to the identity G* = G AQ (1 + 7 y A0 /2), 
which derives from the definitions fl2.21a| ) and ( |2.21b| ) . We have thus recovered without any 



explicit calculation that the deflection of light and its time-delay are both proportional to 
2 + 7 A0 at the 1PK level, and to 2 + 7 at the 1PN level. [Remember that 7 is usually denoted 
as 7 — 1 in the literature, so that 2 + 7 is the usual factor I + 7.] The corrections appearing at 
order 0(Gm/rc 2 ) 2 are due to the five remaining diagrams of Fig. O. Three of them do not 
involve any scalar field, and are thus proportional to G 2 = G ao Gbo(^ + X4 /2)(l + 7 B0 /2). 
The other two involve a scalar line between two material bodies, and yield therefore con- 
tributions proportional to G 2 (aAOtB)$ = — Glj AB /(2 + 7 f AB ), where G\ should be rewritten 
as above in terms of effective gravitational constants. In the case of a single material body 



A, it should be noted that the first two V diagrams of Fig. |TT] do not contribute, since 
the same body cannot be directly connected by a propagator. On the contrary, the first T 
diagram of this figure does contribute, because the two blobs representing the same body A 
are connected to a scalar-scalar-graviton vertex, and not directly to each other. Hence this 
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diagram yields a contribution proportional to 7 AA j (2 + j AA ) in second-order light-deflection 
and time-delay experiments. In conclusion, our diagrammatic approach has allowed us to 

A 

prove is a streamlined way that these experiments do not depend on (3 BC at the 2PK level, 
and more precisely that they can differ from general relativity only by terms of the form 
indicated in Eq. ( |(J1| ) in the case of a single material body A. 

Let us now give explicit expressions for the 2PN light deflection and perihelion advance. 
Following Weinberg pE^] , we use here Schwarzschild-like coordinates {i.e., an "area radius" 
r) and write the integral giving the polar angle in the plane of the trajectory as 



± 



JVABdr/r 2 



[S 2 /c 2 - (m 2 c 2 + J 2 /r 2 )B]2 



(C3) 



where S denotes the conserved energy (including the rest mass contribution) of a test particle 
of mass mo, and J its angular momentum. [Note that these conserved quantities are related 



in a non-obvious way to the constants and Jw used by Weinberg in Eq. (8.4.30) of J49 
_E W = (moc 2 /£) 2 ; Jw = Jc/£] The area radius r should not be confused with the Just 
radius nor the g*-harmonic radius, both denoted also by r in subsection IV-B. To rewrite 
the metric ( |4 . 1 8| ) in Schwarzschild coordinates 



ds 2 /A 2 (ip ) = -B(r)c 2 dt 2 + A(r)dr 2 + r 2 (d6 2 + sin' 



(C4) 



we need to express the area radius r in terms of the isotropic one p. We find easily 

2 



r = p 



i + £±(i + i M ) + -(£±) (1 + 3 r - 47^ - 2? MJ 

per 4 \pc z 
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(C5) 



where p A = G A om A is the Keplerian mass of the attracting body A, and 7 A0 , T are the 
deviations from general relativity in the spatial isotropic metric ( |4.18|) . Equations (|4.19|) 
and (|4.20|) give the expressions of these coefficients in tensor-multi-scalar theories. Note, 
however, that our present calculation is valid for any static and spherically symmetric met- 
ric of the form ( 4.18 ), not necessarily the one predicted by tensor-scalar theories. The 
replacement of ( |C5| ) in (|4.18|) yields 
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The polar angle (|C3|) can now be obtained by a straightforward integration. The case of 
light corresponds to m = 0, and the deflection angle A0 is found to have the form 



= A0! + A0 2 + 0(l/c 6 ) 



(G7) 



where 
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A0! 



2/m(2 + 7ao) 
Poc 2 



//a(2 + 7ao) 



Poc z 



(C8) 



is the 1PK result up to a global correcting factor, and 



7T 

4 



Poc 2 



15 + 3r - 4/3^ + 87^ 



(C9) 



is the actual 2PK contribution. Note that, after having performed the integral (|C3| ) in 
Schwarzschild coordinates, we have expressed the results in terms of po, the minimal dis- 
tance between the light ray and the center of body A, measured in isotropic coordinates. 
[Note that po differs from the impact parameter.] The results are unchanged if one uses 
g*-harmonic coordinate to measure this minimal distance, since the transformation ( |4.22| ) 
introduces corrections only at order 0(l/c 6 ). By contrast, the use of Schwarzschild coordi- 
nates transforms ( |C5| ) into A0 X = [2/i^(2 + 7~Ao)/ r o c2 ] x [1 — (pA/ r o c2 )L while (OJ) remains 
unchanged. This second-order light deflection (p7|)-(|C9|) agrees with previous calculations 



in the literature [0,0,|TH]. It can now be particularized to the case of tensor-multi-scalar 



o 



theories by using the expression (|4.19b|) for T. We find that the term f/3 AA in T cancels 



exactly the —Af3 AA contribution in 



and we get 



A0 2 
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(CIO) 



Together with Eq. ( |C8[ ) , this result confirms therefore the conclusion ( |CT| ) of our diagram- 
matic analysis: The first order deviation from general relativity is proportional to 7^0, while 
the second order contributions involve both 7 A0 and J AA , but not (3 AA . Note also the ap- 
pearance of the same coefficient /ia(2 + 7ao) = 2G>m J 4 in ( |C8| ) and (|C10|) , which is in fact 
a mere rewriting of the bare gravitational constant in terms of observable quantities. 
The coefficient ^ AA /{2 + j AA ) = -(0^)0 entering ( |C10| ) is due to the first T diagram of 



Fig. [IT], where both blobs represent body A. 

Let us now take the second post-Newtonian limit of the above 2PK results, as appropriate 
for interpreting future high-precision experiments in the solar system. The coefficient 7 A0 
of the leading term in ( |C8| ) should therefore be expanded as in Eq. (|5.6|), together with the 
renormalization (|B3b|) due to the gravitational potential of external masses, while the 2PK 
contribution (03) becomes 
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(Cll) 



The explicit expression of the second-order time-delay has been derived in [|T(J for a 
general metric of the form ( |4.18| ). This result confirms also the diagrammatic analysis 
discussed at the beginning of the present appendix: The parameter (3 AA appears again in 
the combination 3r — 4:/3 AA , like in ( |C9|) above, and it vanishes therefore in the case of 
tensor-multi-scalar theories, for which T 



The integral (|C3|) can also be used to derive the periastron advance of a test mass 
m m A . Let us introduce as in Ref. M the notation E = (£ — m c 2 )/m for the specific 
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conserved energy minus rest mass, and h = jT/mofiA for the reduced conserved angular 
momentum. A straightforward integration yields the 2PK periastron shift per orbit: 
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(C12) 



This expression agrees with the general relativistic result derived in [11,9|. Note that con- 
trary to the light-deflection and time-delay formulae, the periastron advance involves not 



only (3 AA but also the 2PK parameter e AAA , entering B = fe^^ + 0([3,j), cf. Eq. ( |4.19a| ). 
The 2PN limit of Eq. ( |C12|) can be obtained easily by using the expansions ( 4.23 ) of the 
different body- dependent parameters. 

The result (|C12| ) is coordinate-independent since it is expressed in terms of the conserved 
quantities E and h. It can nevertheless be helpful to rewrite it in a particular coordinate 
system. Both isotropic and g*-harmonic coordinates give the same result at the 2PK order. 
Let us denote by a(l — e) the coordinate periastron radius, and by a(l + e) the coordinate 
apoastron radius (in isotropic or, equivalently, g*-harmonic coordinates). The conserved 
quantities can then be rewritten as 
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We recover in particular the standard 1PN formula, A0 = 2irfi A (3 — (3 + 27)/a(l — e 5 



)c 2 + 



0(l/c 4 ), in terms of the semilatus rectum a( 
in Eq. (|5.9| ) above. 



1 — e ). The 2PN result takes the form displayed 
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FIGURES 

FIG. 1. Diagrammatic interpretation of the effective gravitational constant G = G*Aq(1 + «g). 

FIG. 2. Diagrammatic interpretation of the contraction (a/5a)o involved in the Eddington 
parameter (3 = ^(a/3a)o/(l + «o) 2 - 

FIG. 3. Diagrammatic notations for the material sources, the fields, and their propagator. 

FIG. 4. Diagrammatic expression of the 3> l -linear terms of the total action ( |3.6| ), for i = 1,2,3,4. 

FIG. 5. Equation ( jTgg ) satisfied by the field ¥[er]. 



FIG. 6. Diagrammatic expansion of the Fokker action ( |3.13D . 

FIG. 7. Expression of the diagrams of Fig. |6| when the graviton and scalar propagators are 
represented respectively as curly and straight lines. 

FIG. 8. Constraints imposed by four different binary-pulsar data on the 2PN parameters e, C 

FIG. 9. Region of the e— £ plane allowed at the la level by the four tests of Fig. |8|. 

FIG. 10. Decomposition of the F and H diagrams of Fig. |(| showing in bold the lines involving 
derivatives in the vertex <&d<&d&. Such a decomposition is useful to compute the last two F and H 
diagrams of Fig. [?] without symmetrizing the 3-graviton vertex they contain. 

FIG. 11. Diagrams describing the interaction of light (wavy lines) with the gravitational field 
generated by material bodies, at orders 0(Gm/rc 2 ) -1PK- and 0(Gm/rc 2 ) 2 -2PK-. 
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TABLES 



TABLE I. Expression of Nordtvedt's mass parameters in tensor-multi-scalar theories of gravity. 

Nordtvedt's parameters Tensor-multi-scalar theories 

WT) f^A 

T ij /M{I) i M{I) j G AB 

7©ij/Lij 1 + 

(2/3 - \)Y ijk IM{I)iM{I)jM{I) k (1 + 2Jbc)GabGac 

M{G)/M(I) G A0 /G 

M( 7 )/M(J) G A0 (1 + Tao)M1 + 7)_ 

M(J3)/M(I) G A0 (1 + 2P A0 )/G(1 + 2/3) 

M((3')/M(I) G\ (l + 2^ )/G 2 (l + 2£) 
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Figure 10 
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